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Resume en frangais 



Dans cette these, nous etudions les proprietes statiques et dynamiques des echelles 
de spin-1/2 soumises a un champ magnetique. Faiblement couplees, ces echelles 
permettent d'etudier a la fois la physique du liquide de Luttinger (LL) apparais- 
sant dans de nombreux systemes unidimensionnels (ID) et la condensation de Bose- 
Einstein (BEC) qui est un effet typiquement tridimensionnel (3D). Notre travail a ete 
en grande partie motive par le compose (C5Hi2N)2CuBr4 (BPCB) recemment syn- 
thetise. Ce materiau est considere conmie ay ant une simple structure d' echelles cou- 
plees oil des interactions plus complexes (frustration, anisotropic ou Dzyloshinskii- 
Moriya) restent faibles. De plus, ses couplages d'echange sont suffisamment faibles 
pour rendre 1' ensemble de son diagranmie de phase experimentalement accessible 
en appliquant un champ magnetique. 

Pour etudier ces systemes, nous utilisons une combinaison de methodes analy- 
tiques (theorie du liquide de Luttinger et technique de bosonisation) et numeriques 
(groupe de renormalisation de la matrice densite (DMRG)). La premiere est une 
theorie des champs permettant de decrire la physique de basse energie de nombreux 
systemes ID sans bande interdite telle que la phase aimantee des echelles de spin- 
1/2. La seconde est une methode variationnelle particulierement bien adaptee aux 
systemes ID et permettant de calculer leurs proprietes a temperature nulle. Elle per- 
met egalement d'extraire les parametres du LL a partir du calcul des correlations sta- 
tiques et de I'aimantation pour obtenir une description quantitative de la physique de 
basse energie. La methode DMRG a ete recemment etendue au calcul des proprietes 
a temperature finie ainsi que leur evolution temporelle. Cette demiere extension 
est notamment utilisee, dans ce travail, pour calculer les fonctions de correlation 
dynamiques. Le faible couplage inter-echelles est, quant a lui, pris en compte en 
utilisant une approximation de champ moyen. 

Dans un premier temps, nous explorons le diagramme de phase des echelles non- 
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couplees en etudiant leurs proprietes thermodynamiques. Nous calculons leur mag- 
netisation et leur chaleur specifique. Ces deux quantites revelent des caracteristiques 
de basse temperature fideles a la description du LL. Au contraire, leur comportement 
a haute temperature est principalement dicte par les excitations de triplets a haute en- 
ergie. On determine egalement a partir des extrema de ces deux quantites la limite 
de validite approximative de la description du LL. La comparaison de nos calculs 
avec les mesures effectuees sur BPCB dans un domaine complet de temperature et 
champ magnetique est excellente. Ceci confirme la simple structure d'echelle de ce 
compose. De plus, le domaine de validite de la description du LL pour BPCB est 
un ordre de grandeur plus grand que la temperature de transition de I'ordre 3D a 
basse temperature. Ce qui laisse ainsi un large domaine de temperature pour tester 
le liquide de Luttinger. 

Pour ce faire, nous calculons les predictions du LL du temps de relaxation mesure 
par resonance magnetique nucleaire (NMR). De plus, la prise en compte du couplage 
inter-echelles par une approximation de champ moyen permet d'acceder, a I'aide 
de la description du LL des echelles isolees, a la temperature critique et I'ordre 3D 
transverse antiferromagnetique associe. Ces demiers caracterisent la transition BEC. 
La mesure de ces trois quantites a I'aide d' experiences NMR et de diffraction de 
neutrons sont toutes en tres bon accord avec les predictions du LL. EUes permettent 
done de tester trois differentes fonctions de correlations calculees a partir des memes 
parametres du LL. Ceci foumit le prenuer test quantitatif de la theorie du liquide de 
Luttinger. 

Dans un second temps, nous avons etudie les fonctions de correlation dynamiques 
a temperature nuUe des echelles non couplees. Les excitations fournissent d'importantes 
informations sur le systeme et permettent ainsi de le caracteriser en detail. En partic- 
ulier, nous presentons leur interessante evolution avec le champ magnetique applique 
et pour differents couplages. Le continuum apparaissant a basse energie est qualita- 
tivement decrit a I'aide d'une approximation par une chaine de spin- 1/2 anisotrope 
equivalente a 1' approximation de fort couplage de I'echelle. Cette approximation 
n'est en revanche pas valable pour la description des excitations de moyenne et haute 
energie. En effet, ces demieres necessitent la prise en compte des triplets de haute en- 
ergie negliges par cette approximation. Fait interessant, les excitations de moyenne 
energie peuvent etre decrites par un modele t-J et presentent done des caracteristiques 
typiques des systemes itinerants. On verifie de plus que revaluation numerique de 
ces correlations valable a moyenne et haute energie convergent correctement sur la 
description donnee par le LL restreint aux excitations de basse energies. 

Les mesures de diffusion de neutron inelastique (INS) etant directement reliees 
aux correlations dynamiques, on fournit une prediction complete des spectres mesures 
sur BPCB. II est gratifiant de noter que la resolution en energie et quantite de mou- 
vement de nos calculs est actuellement meilleure que celle des experiences. Les 
mesures se limitant pour 1' instant aux excitations de basse energie, il est difficile d'y 
distinguer la difference avec la prediction fournie par une echelle de spin- 1/2 et celle 
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donnee par 1' approximation de fort couplage. Des mesures du spectre de moyenne 
et haute energie comportant des excitations caracteristiques du modele sous-jacent 
permettraient de raffiner 1' etude experimentale sur BPCB. 

Plus generalement, d'un point de vue conceptuel ou en lien avec BPCB plusieurs 
points necessitent une etude plus etendue. 

La prise en compte du couplage inter-echelles et de la temperature pour les pre- 
dictions theoriques devrait etre etendue aux quantites dynamiques. En effet, il serait 
interessant d'etudier leur impact sur les excitations du systeme. Des phenomenes 
tels que le deplacement ou I'elargissement des excitations avec la temperature ont ete 
observes dans d'autres systemes magnetiques. Une etude detaillee pres des champs 
critiques serait particulierement recommandee. Dans ces regimes, le systeme subit 
une transition dimensionnelle entre un etat ID (ou les excitations sont essentielle- 
ment fermioniques) et un etat 3D (ou les excitations ont une description bosonique). 
De plus, il serait important de prendre en compte la structure tridimensionnelle reelle 
du compose. Ceci nous permettrait egalement de comprendre plus en detail les devi- 
ations de I'ordre 3D mesure experimentalement sur BPCB avec notre prediction util- 
isant une approximation de champ moyen. 

Recemment, de faibles anisotropics ont ete detectees sur BPCB a I'aide de mesures 
de resonance electron- spin. Les effets de ces anisotropics sur la physique des echelles 
etant actuellement peu connue, une etude approfondie de ces phenomenes serait per- 
tinente. Elle permettrait notanmient de comprendre certaines deviations entre les 
mesures sur BPCB et nos predictions theoriques utilisant des echelles isotropes. 

Recemment des singularites apparaissant dans les correlations dynamiques et 
sortant de la description du LL ont ete mises en evidence. Or, pour qu'une etude 
de ces effets sur les excitations des echelles soit possible, il faudrait ameliorer la 
precision des correlations calculees numeriquement. L' optimisation des algorithmes 
existant ou le developpement de nouvelles methodes plus performantes pour attein- 
dre ce but fait done partie integrante des extensions futures. 

Etant donne que le compose BPCB est maintenant bien modelise par de sim- 
ples echelles de spin- 1/2, il serait interessant de lui aj outer des impuretes pour y 
etudier des phenomenes plus complexes. En effet, la presence de desordre perme- 
ttrait r etude du verre de Bose alors que I'ajout de porteurs de charges pourraient 
faire apparaitre un etat supraconducteur. D'un autre cote, notre demarche etant assez 
generale, elle pourrait etre etendue a 1' etude d'autres composes d'echelle tels que 
DIMPY recemment synthetise ou des structures plus complexes. 
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Abstract 



We investigate weakly coupled spin-1/ 2 ladders in a magnetic field. The work is mo- 
tivated by recent experiments on the compound (C5Hi2N)2CuBr4 (BPCB). We use 
a combination of numerical and analytical methods, in particular the density matrix 
renormalization group (DMRG) technique, to explore the phase diagram and the ex- 
citation spectra of such a system. We give detailed results on the temperature depen- 
dence of the magnetization and the specific heat, and the magnetic field dependence 
of the nuclear magnetic resonance (NMR) relaxation rate of single ladders. For cou- 
pled ladders, treating the weak interladder coupling within a mean field approach, we 
compute the transition temperature of triplet condensation and its corresponding an- 
tiferromagnetic order parameter. Existing experimental measurements are discussed 
and compared to our theoretical results. Furthermore we compute, using time depen- 
dent DMRG, the dynamical correlations of a single spin ladder. Our results allow 
to directly describe the inelastic neutron scattering cross section up to high energies. 
We focus on the evolution of the spectra with the magnetic field and compare their 
behavior for different couplings. The characteristic features of the spectra are inter- 
preted using different analytical approaches such as the mapping onto a spin chain, 
a Luttinger liquid (LL) or onto a t-J model. For values of parameters for which such 
measurements exist, we compare our results to inelastic neutron scattering experi- 
ments on the compound BPCB and find excellent agreement. We make additional 
predictions for the high energy part of the spectrum that are potentially testable in 
future experiments. 
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CHAPTER 1 



Introduction 



In many condensed matter systems the quantum fluctuations and the interactions 
between particles play a crucial role. Various important effects such as the high 
temperature superconductivity, the fractional quantum Hall effect or Mott insulators 
arise in the so-called strongly correlated quantum systems. In these systems the in- 
teraction is comparable to other energy scales and thus needs to be treated on equal 
footing. In a Mott insulator, due to the Pauli principle, the interplay between interac- 
tions and kinetic energy can induce a strong antiferromagnetic spin superexchange. 
Such exchange leads to a remarkable dynamics for the spin degrees of freedom. On a 
simple square lattice, the antiferromagnetic exchange can stabilize an antiferromag- 
netic order. By variations in dimensionality and connectivity of the lattice a variety 
of complex phenomena can arise |(T|, for instance, spin liquid ||2|, Bose-Einstein 
condensation |j3||5J (BEC), Luttinger liquid (6^ J] (LL) or Haldane gap |8|. Recently, 
among those effects two fascinating situations in which the interaction strongly fa- 
vors the formation of dimers have been explored in detail. 

The first situation concerns a high dimensional system in which the antiferromag- 
netic coupling can lead to a spin liquid state made of singlets along the dimers. In 
such a spin liquid the application of a magnetic field leads to the creation of triplons 
which are spin-1 excitations. The triplons which behave essentially like itinerant 
bosons can condense leading to a quantum phase transition that is in the universality 
class of BEC. Such transitions have been explored experimentally and theoretically 
in a large variety of materials, belonging to different structures and dimensionali- 
ties f9\ . On the other hand, low dimensional systems behave quite differently. Quan- 
tum fluctuations are extreme, and no ordered state is usually possible. In many quasi 
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one-dimensional systems the ground state properties are described by LL physics 
that predicts a quasi long range order. The elementary excitations are spin- 1/2 ex- 
citations (spinons). They behave essentially as interacting spinless fermions. This 
typical behavior can be observed in spin ladder systems in the presence of a mag- 
netic field. Although such systems have been studied theoretically intensively for 
many years in both zero 1 10 - 18] and finite magnetic field |[3|[T9}|29j, a quantitive 
description of the LL low energy physics remained to be performed specially for a 
direct comparison with experiments. 

Quite recently the remarkable ladder compound fSOl (C5Hi2N)2CuBr4, usually 
called BPCB (also known as (Hpip)2CuBr4), has been investigated. The compound 
BPCB has been identified to be a very good realization of weakly coupled spin lad- 
ders. The fact that the interladder coupling is much smaller than the intraladder 
coupling leads to a clear separation of energy scales. Due to this separation the 
compound offers the exciting possibility to study both the phase with Luttinger liq- 
uid properties typical for low dimensional systems and the BEC condensed phase 
typical for high dimensions. Additionally, the magnetic field required for the re- 
alization of different phases lies for this compound in the experimentally reachable 
range. Actually various experimental techniques such as nuclear magnetic resonance 
(NMR) [31^1, neutron diffractior[^ (ND) |32|, specific heat and magnetocalorific ef- 
fect p3| are used to probe the static properties of different phases of this compound. 
In order to interpret correctly these experiments a quantitative theoretical description 
of weakly coupled spin ladders is thus strongly required. 

The excitations of this compound have recently been observed by inelastic neu- 
tron scattering | 34'351 experiments (INS). These are directly related to the dynamical 
correlations of spin ladders. Although these dynamical correlations have been inves- 
tigated intensively in absence of a magnetic field during the last decades p4}|T8p6| , 



a detailed analysis and a quantitative description of their magnetic field dependence 
specially for the high energy excitations is clearly missing. The direct investigation 
of such excitations is of high interest, since they not only characterize well the spin 
system, but the properties of the triplon/spinon excitations are also closely related to 
the properties of some itinerant bosonic/fermionic systems. Indeed using such map- 
pings |[6j of spin systems to itinerant fermionic or bosonic systems, the quantum spin 
systems can be used as quantum simulators to address some of the issues of itinerant 
quantum systems. One of their advantage compared to regular itinerant systems is 
the fact that the Hamiltonian of a spin system is in general well characterized, since 
the spin exchange constants can be directly measured. The exchange between the 
spins would correspond to short range interactions, leading to very good realization 
of some of the models of itinerant particles, for which the short range of the interac- 
tion is usually only an approximation. In that respect quantum spin systems play a 
role similar to the one of cold atomic gases |37|, in connection with the question of 



^ND consists in elastic neutron scattering by opposition INS implies an energy transfer. This technique 
can be used to measure the long range magnetic orders and is shortly discussed in Sec. |5. 6.1] 
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itinerant interacting systems. 

In this thesis, we present an analysis of the properties of weakly coupled spin- 
1/2 ladders in a magnetic field. We consider both the low energy physics and the 
excitations providing a quantitative description necessary for an unbiased compari- 
son with experiments. The main achieved results discussed in this thesis as well as 
in Ref. ||3T}f3|[38| are: 

• Combining a LL analytical technique and numerical density matrix renormal- 
ization group (DMRG) methods, we provide a quantitative description of the 
static and dynamic properties of spin- 1/2 ladders in a full range of temperature 
and energy. 

• We provide a detailed analysis of the dynamical correlations in a full range of 
magnetic field and couplings. 

• Taking into account a weak interladder coupling by a mean field approxima- 
tion, we characterize the BEC of triplons appearing at low temperature in 
weakly coupled spin- 1/2 ladders. 

• Comparing the experimental measurements on the compound BPCB to our 
theoretical computations, we confirm the weakly coupled spin-1 / 2 ladder struc- 
ture of this material which provides the first quantitative test of the LL theory 
and shows a phase transition to a BEC of triplons. 



1.1 Plan of the thesis 



Although strongly based on Ref. | |38| , this thesis contains various important technical 
aspects as well as introductions on the methods and broader discussions which were 
omitted in Ref. 1 3^. In addition it includes several comparisons with the experiments 
on BPCB from Ref. (3TH33J. The plan of the thesis is as follows. 



In chapter [2j we define the model of weakly coupled spin ladders. Its basic 
excitations and phase diagram are introduced as well as the spin chain map- 
ping which proves to be very helpful for the physical interpretations. The spin 
ladder compound BPCB is also characterized with a detailed discussion of its 
chemical structure as well as the resulting interactions. The chapter [2]is a gen- 
eral introduction on spin- 1/2 ladder systems and its experimental realizations. 
It can be easily skipped by an informed reader. 

The chapter [3] provides a description of various theoretical techniques suited 
for dealing with low dimensional systems. We introduce the DMRG methods 
as well as the LL theory focussing on their application on spin- 1/2 ladder. 
In particular, we introduce the recent real-time variant of DMRG to obtain 
the dynamics |39 -42| in real time and the dynamical correlation functions. 
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A similar technique is also presented to obtain finite temperature results (43]- 
[45| . The effect of a weak interladder coupling is discussed using a mean field 
approximation. The chapter [3] is technically oriented and presents how the 
achieved results shown in |4] and [5] are computed. 

• In chapter |4] we give a detailed characterization of the phase diagram of 
weakly coupled spin- 1/2 ladders focusing on their static properties (magneti- 
zation, specific heat, BEC critical temperature, order parameter) and the NMR 
relaxation rate. This characterization is followed by a comparison with the 
measurements on BPCB. 

• The chapter |5] presents the computed dynamical correlations of a single spin 
ladder at different magnetic fields and couplings. The numerical calculations 
are compared to previous results (linked cluster expansion, spin chain map- 
ping, weak coupling approach) and analytical descriptions (LL, t-J model). 
The theoretical spectra are compared to the low energy INS measurements on 
the compound BPCB and provide predictions for the high energy part of the 
INS cross section. The effects of a low interladder coupling on the dynamics 
are briefly discussed as well as the ND technique for measuring long range 
magnetic order. 

• In chapter|6j we summarize our results and discuss further perspectives. 

• The appendix [Ajpresents the strong coupling expansion which provides a sim- 
plified picture for the understanding of many results presented in this thesis. 

• The appendix |b] summarizes the static correlations in spin- 1/2 chains and lad- 
ders computed in Ref. | [28| for finite size systems. These enter as a key point 
in the quantitative description of the low energy physics of spin- 1/2 ladders. 
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Spin-1/2 ladders 



The physics of quantum spin systems depends strongly on their microscopic char- 
acteristics. In fact, a large variety of phenomena can appear depending on the local 
spin S, the type of interaction and the geometry of the system. The external envi- 
ronment is also very important. For instance, applying a magnetic field or pressure 
on the system can lead to even richer physics. Various combinations of internal con- 
straints as well as external conditions have been investigated theoretically for many 
years These studies have lead to several fundamental discoveries such as the 
specific properties of ID spin chains which are gapped or gapless in case the local 
spins S is integer or half-integer 1 8 1 respectively or the strong connection between 
the high temperature superconductors and the 2D quantum magnetism on a square 
lattice (461. 

Motivated by these two important results, the spin-1/2 ladders lying between 
these two ID and 2D limits, have been studied theoretically intensively in both 
zero (TOjITSj and finite magnetic field |^[T9jj29|. In addition the identification of 
the (C5Hi2N)2CuBr4 (BPCB) compound [30T as a very good realization of weakly 
coupled spin-1/2 ladders has even increased the interest in these systems. Indeed 
due to its particularly low energy couplings its complex phase diagram (shown in 
Fig. |2.2[ b) is fully accessible experimentally by tuning a magnetic field and has been 
explored with various techniques [30-^35^^38^47 - 50]. 

In this chapter we first describe the weakly coupled spin-1/2 ladder model on 
which we will focus in this work. Next we remind briefly its main physical features 
and introduce the spin chain mapping which provides a simple interpretation for 
certain features. Finally we present the compound BPCB with an analysis of its 
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chemical structure and the resulting interactions. 



2.1 Weakly coupled spin- 1/2 ladders 

Various spin- 1/2 ladder systems have been investigated including different coupling 
geometries with frustration or long range interaction as well as site dependent or 
anisotropic interactions (12 19 27 5T|. In this work we consider a simple ladder 



structure with isotropic Heisenberg couplings between nearest neighbors and no frus- 



2.1 



tration as pictured in Fig. 2.1 a. In addition, a weak interladder coupling is discussed 
with the assumed unfrustrated 3D coupling structure shown in Fig. |2.1[ b. The gen- 
eral Hamiltonian for these weakly coupled spin- 1/2 ladders is 

H3D = ^ Hjj^ + J' ^ Si^k,fi ' ^u,k',fi' (2.1) 

where is the Hamiltonian of the single ladder /i and J' is the strength of the 
interladder coupling. The operator Si^k,n ~ i'^f,k,fi^ '^ik /i' ^i,k,^) ^^^^ ^^^^ ^ 

(/ = 1, 2, . . . , L) of the leg /c (/c = 1, 2) of the ladder /i. Often we will omit ladder 
indices from the subscripts of the operators (in particular, replace Si^k,^ with Si^k) 
to lighten notation. S^^ (a = are conventional spin- 1/2 operators with (we 

mostly use H = ks = I) 

[Sr,k,S^^,}=ie^MSl„ and Sf;, = Sf^, ± iSl, (2.2) 

P^j = x^y^ z and ea(3-f is the totally antisymmetric tensor. 
The spin-1 / 2 Hamiltonian of the spin-1 / 2 two-leg ladder illustrated in Fig 

is 

H^ = J^H^^J\\H\\ (2.3) 
where J±_ (Jy) is the coupling constant along the rungs (legs) and 

H^ = Y. • - h'J^'M^ H\\ = S,,, • S,+i,, (2.4) 

I l,k 

The magnetic field, , is applied in the z direction, and is the 2:-component of 
the total spin operator M = (S/,i+S/,2). Since has the symmetry —h^, 
— M^, we only consider > 0. The relation between and the physical 
magnetic field in experimental units is given in Eq. ( |2.16| ). 

In this work we focus on the case of spin-1 /2 antiferromagnetic ladders weakly 
coupled to one another. This means that the interladder coupling J' > is much 
smaller than the intraladder couplings Jy and J±, i.e. 

< J' < J|| and J^. (2.5) 

Therefore, the interladder coupling J' will be treated perturbatively by a mean field 
approximation (see Sec. |3.3| ) neglecting the microscopic details of the interladder 
interactions pictured in Fig.|2.1[b for the supposed coupling structure of BPCB. 
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(a) !✓ !> (b) l> V W ^ 

'7,' ! : ' 7/ ' i ' 

^± W I ; I J I i I 

II , : I I ; I 

II , i I I ; I 

II ,111 

Figure 2.1: (a) Schematic representation of a single ladder: 8^,^ are the spin operators 
acting on the site / of the leg /c = 1, 2. (b) Schematic representation of weakly coupled 
ladders (with the supposed interaction geometry of BPCB shown in Fig. |2.6| l: Si,fc,/x is 
the spin operator acting on the site / of the leg /c = 1, 2 of the ladder \i. In (b), the 
ladders are oriented perpendicular to the sheet. J±, J|| and J' are the coupling along the 
rungs and the legs, and the interladder coupling represented by thick, medium and thin 
dashed lines, respectively. The dotted rectangle represents an unit cell. 



2.1.1 Spin ladder to spin chain mapping 



The physical properties of a single ladder ( |2.3| ) are defined by the value of the di- 
mensionless coupling 

^11 

7=/. (2.6) 

In the limit Jy = (therefore 7 = 0) the rungs of the ladder are decoupled. We 
denote this decoupled bond limit (DBL) hereafter. The four eigenstates of each de- 
coupled rung are: the singlet state 



in) - lit) 



(2.7) 



with the energy Es = — 3J^/4, spin 5 = 0, and z-projection of the spin = 0, 
and three triplet states 

|t+) = in), 10 = ^^^^^^, r) = iu) (2.8) 

with S' = 1, 6'^ = 1, 0, -1, and energies E^^ = - h\ Eto = Et- = 

Jx/4 + h^, respectively. The ground state is \s) below the critical value of the 
magnetic field, h^^^ = J±, and above. The dependence of the energies on the 
magnetic field is shown in Fig.|2.2[a. 



A small but finite 7 > delocalizes triplets and creates bands of excitations with 
a bandwidth ^ Jy for each triplet branch. This leads to three distinct phases in the 
ladder system p.3|) depending on the magnetic field: 
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Figure 2.2: (a) Energy of the triplets |t^), |t°), \t~) (solid lines) and singlet \s) (dashed 
line) versus the applied magnetic field in the absence of an interrung coupling (J\\ =0). 
The dotted lines represent the limits of the triplets excitation band when J\\ ^ 0. (b) 
Phase diagram of weakly coupled spin ladders: crossovers (dotted lines) and phase tran- 
sition (solid line) that only exists in the presence of an interladder coupling are sketched. 
[Taken from Ref. |38|] 



(i) Spin liquid phas^ which is characterized by a spin-singlet ground state (see 



Sec. |4.1[ ) and a gapped excitation spectrum (see Sec. |5.2| ). This phase appears 
for magnetic fields ranging from to hc^ . 

(ii) Gapless phase, which is characterized by a gapless excitation spectrum. It 
occurs between the critical fields hc^ and hc^ . The ground state magnetization 
per rung, = {M^) /L, increases from to 1 for running from hc^ to 



/ic2 • The low energy physics can be described by the LL theory (see Sec. 3.2 ). 



(iii) Fully polarized phase, which is characterized by the fully polarized ground 
state and a gapped excitation spectrum. This phase appears above hc2 • 

The transition between (i) and (ii) can also occur in several other gapped systems 
such as Haldane 5* = 1 chains or frustrated chains | [2T||52}|54| . In the gapless phase, 



^This phase is also called quantum disordered |9 1. It appears in other antiferromagnetic systems such 
as frustrated antiferromagnets. 
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the distance between the ground state and the bands \t^) and which is of the 
order of Jj_, is much larger than the width of the band ~ Jy , since 7 <C 1. 

For small 7 the ladder problem can be reduced to a simpler spin chain problem. 
The essence of the spin chain mapping |3 , 19, 55 , 56J is to project out \t^) and \t~) 
bands from the Hilbert space of the model ^23]^ The remaining states \s) and 
are identified with the spin states 

\l) = \s), |f) = |i+). (2.9) 

The local spin operators S^^^ can therefore be identified in the reduced Hilbert space 
spanned by the states ( |2.9| ) with the new effective spin-1/2 operators S^: 

^<^ = ^^^ 5;% = i(l + 2^f). (2.10) 

The Hamiltonian ( |2.3| ) reduces to the Hamiltonian of the spin-1/2 XXZ Heisenberg 
chain 

Hxxz = J|| E (^f^f+i + ^f^li + AS!SU,)-h^M^^L (-^ + y - y) • 

(2.11) 

Here the pseudo spin magnetization is — SI, the magnetic field = — 
-^_L — -^||/2 and the anisotropy parameter 

A = i. (2.12) 

Note that the spin chain mapping constitutes a part of a more general strong coupling 
expansion of the model ( |2.3| ), as discussed in the appendix [A| 
For the compound BPCB the parameter 7 is rather small 

7«^« 0.282. (2.13) 

and the spin chain mapping \2Al\ gives the values of many observables reasonably 
well. Some important effects in particular at high energy are, however, not captured 
by this approximation. For instance, due to their connection with high energy triplet 
excitations, several correlations cannot be described by this approximation. Other 
examples will be given in later chapters. 

2.1.2 Role of weak interladder coupling 

Let us now turn back to the more general Hamiltonian \2A) and discuss the role of 
a weak interladder coupling J' . The spin liquid and fully polarized phase are al- 
most unaffected by the presence of J' whenever the gap in the excitation spectrum 
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is larger than J' (see, e.g., Ref. pO| for more details). However, a new 3D anti- 
ferromagnetic order in the plane perpendicular to emerges in the gapless phase 
for T ^ J^ The corresponding phase, called SD-ordered, shows up at low enough 
temperatures Tc in numerous experimental systems with reduced dimensionality and 
a gapless spectrum ||9|. This phase transition is in the universality class of Bose- 
Einstein condensation and is discussed in more detail in Sees. 13.31 and 14.41 For the 
temperature T > J' the ladders decouple from each other and the system undergoes 
a deconfinement transition into a Luttinger liquid regime (which will be described in 
Sec.|3.2|). For T > Ju the rungs decouple from each other and the system becomes 



a (quantum critical) paramagnet. The transition from the 3D ordered to the LL [24 | 
and the crossover from the LL to the quantum critical regime | [23| induce specific 
features in several thermodynamic quantities such as the specific heat or the mag- 
netization. These characteristics are pointed out in Sees. |4.2.T] and [4Z!2| and used 
to locate the LL to quantum critical crossover. All the above mentioned phases are 
illustrated in Fig.|2.2[b. 



2.2 Experimental realizations of spin-1/ 2 ladders 



The spin- 1/2 ladder structure shown in Fig. |2.1[ a has been pointed out in several 
materials containing Cu^^ ions with an unpaired external electronic orbital. Ini- 
tially motivated by their strong connection with high temperature superconductors 
the inorganic compounds such as SrCu2 03 |57 | or La4SrioCu2404i | 58 | have been 
intensively investigated during the 90' s. Although these materials show typical fea- 
tures of spin- 1/2 ladders, their strong antiferromagnetic superexchange couplings 
(^ 1300 K) occurring through Cu — O — Cu bonds induces a large spin gap, hd. 
This big spin gap is of the order of few hundreds of Kelvin and thus prevents any 
investigation of the magnetic field effects. 

During the last decade new organic compounds have shown similar spin lad- 
der structures. Mediated by long organic chains, the antiferromagnetic superex- 
change couplings are usually much smaller than these observed in inorganic materi- 
als. Thereby, applying a magnetic field allows one in principle to explore the whole 



phase diagram shown in Fig. |2.2| b which was totally inaccessible experimentally 
for the inorganic materials. First investigated, the compound Cu2(C5Hi2N2)2Cl4 
|[59^!6Q | has finally shown significant deviations from the simple spin ladder structure 



(Fig. |2.1[ a). The presence of frustration or more complicated coupling paths | |6T|[62| 
as well as Dzyaloshinskii-Moriya interactions |[63]|64j is actually debated. 

More recently the compound (C5Hi2N)2CuBr4 first presented in Ref. |30| and 
commonly called BPCB or (Hpip)2CuBr4 has been intensively investigated using 
different experimental methods such as nuclear magnetic resonance | [3T| (NMR), 
neutron diffraction |32| (ND), inelastic neutron scattering ||34||35J (INS), calorime- 
try [33 1, magnetometry (471, magnetostriction ^48^|49J, and electron spin resonance 
spectroscopy jSOj (ESR). Except for small coupling anisotropics [50J , which are 
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briefly discussed in Sees. 4.5 and 2.2 no significant deviations from the simple lad- 
der structure (Fig. |2.1[ a) have been detected. In addition, a small interladder coupling 
J' <C J±, J|| has been pointed out in Refs. | [3H[32|. Although the exact 3D interlad- 
der coupling structure 1 32 1 is discussed in Sec. |4.5[ this interladder coupling allows 
us to explore the transition from the ID to the 3D regime which is experimentally 
accessible in BPCB. This compound is thus an extraordinary experimental tool for 
exploring the weakly coupled spin- 1/2 ladder phase diagram depicted in Fig. |2.2[ 
All the predicted phases have been observed in this compound and Figs. |2.3| and |2.4| 
show the experimental determination of its phase diagram from specific heat and 
neutron diffraction measurements, respectively. A detailed analysis of these exper- 



iments is presented in Sec. 4.5 In this thesis, we focus mainly on this compound 
which provides a strong motivation and an experimental test for our theoretical in- 
vestigation. 

The compound (C5Hi2N)2CuBr4 has the chemical structure represented in Fig. 
2.5 This structure f30l is monoclinic, space group P2i/c f6?l, with the unit cell 



dimensions a = 8.487 A, 6 = 17.225 A and c = 12.380 A (a, b and c are the unit 
ceU vectors of BPCB) and the angle (3 = 99.29° between the a and c axes. 

The magnetic properties of the compound are related to the unpaired highest 
energ y electron ic orbital of the Cu^+ ions. Thus the corresponding spin structure 
(Figs. [2!6|[2!5] and |2l^ matches with the Cu^^ location |[30||3l]|34). The unpaired 



spins- 1/2 interact together by antiferromagnetic superexchange coupling through 
the long organic chains and form two types of equivalent weakly coupled ladders 



(Fig. 2.6) along the a axis. The direction of the rung vectors of these ladders are 

di,2 = (0.3904, ±0.1598, 0.4842) (2.14) 



in the primitive vector coordinates (Fig. 2.6 b). Thus due to their different orientation 



the two types of ladders become slightly distinct when a magnetic field is applied 
breaking the structure symmetry. 

The intraladder couplings from Eq. ( |2.3| ) were determined to be Jx ~ 12.6 — 
13.3 K, J|| ^ 3.3 — 3.8 K with different experimental techniques and at different 
experimental conditions | [30|j35l[47|j49| . In this work, we use the value^ 

J^^12.6K, J||^3.55K. (2.15) 

Recently, a slight anisotropy of the order of 5% of J± has been discovered by 
ESR 1 50 1 measurements. This anisotropy could explain the small discrepancies be- 
tween the couplings found in different experiments. 
The magnetic field in Tesla is related to replacing 

giiBh^ (2.16) 



^These parameters were the first outputs from the NMR measurements, which were later refined to 
the values from Ref. |31 1. Note that small changes in these values do not affect the main results of the 
calculations. 
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4 6 8 10 
Magnetic Field [T] 



Figure 2.3: False color representation of the measured ratio of the specific heat over the 
temperature, c/T, performed on BPCB versus the applied magnetic field, h^, and the 
temperature, T (see also Fig. 4.10 for c(T) plots for several magnetic fields). The ex- 
perimentally determined phase diagram shows the various phases sketched in Fig. |2.2| b: 
spin liquid (called here quantum disorder (QD)), quantum critical (QC), and LL regimes. 
The two critical fields, hd and hc2, are denoted here by Be and Bs. Local maxima from 
the reduction of the triplet gap by the Zeeman effect are indicated by crosses. Cir- 
cles denote the LL crossover based on measurements of the magnetocaloric effect (see 
Fig. [4lT| using the {dm^ /dT) \ 



criterium discussed in Sec. 4.2.1 
blue line indicates the onset of long-ranged order below approximatively 100 mK (3D 
ordered phase in Fig. |2.2| b) shown in more details in Figs. |2.4| and |4.6| In the inset, the 
lattice structure of BPCB in projection along the b axis, is depicted with Cu atoms in 
blue and Br in white. [Taken from Ref. p3)] 



The dashed 



in Eq. \23) with (ib being the Bohr magneton and g being the Lande factor of the 
unpaired copper electron spins. The latter depends on the orientation of the sample 
with respect to the magnetic fielcj^ For the orientation chosen in the NMR measure- 
ments 1 31 1, the Lande factor amounts to g ^ 2.126. It can vary up to ~ 10% for 
other experimental setups pO| . 

As one can see from the projection of the spin structure onto the plane perpen- 



dicular to the a axes (Fig. 2.6 b), each rung is expected to have nc = 4 interladder 



neighboring spins. The interladder coupling J' has been experimentally determined 



^More precisely, the Lande factor is different for each ladder forming the compound and varies with 
their orientation with respect to the magnetic field. 
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Magnetic field [T] 



Figure 2.4: Low temperature phase diagram of BPCB. The crossover temperature to 
the LL regime (black circles) is derived from the magnetocaloric effect represented by 
false color (see also Fig. 4.11 1 using the {dm^/dT)\h^ = criterium discussed in 



Sec. |4.2.l] The phase transition to the 3D ordered phase, Tc, has been detected from the 
appearance of the magnetic order in ND measurements (red circles) (see also Fig. |5.13| l. 
See Fig. |4.6| for a comparison with the NMR determination of Tc. The red line is the 
theoretical prediction of Tc discussed in Sees. 



3.3 



and 



4.4 



[Taken from Ref. ^32J] 




Figure 2.5: Chemical structure of BPCB. The 10 protons attached to the C atoms are 
not shown. Solid thick blue lines and dashed thick blue lines stand for the interaction 
path J± and Jy, respectively. [Taken from Ref. | ^31J ] 



to be |[3Tl[32l| 



f ^ 20 - 100 mK. 



(2.17) 



As we will discuss in Sec. 4.5 the exact 3D coupling structure shown in Fig. 2.6 and 
the precise value of J' are actually debated ||31j^38J . 



13 



2. SPIN-1/2 LADDERS 




Figure 2.6: Coupling structure of BPCB where the unpaired electron spins of the Cu^^ 
atoms in the first (second) type of ladders are pictured by red (blue) spheres. The J^, J\\ 
and J' coupling paths are represented in turquoise, pink, and green, respectively, a, b 
and c are the three unit cell vectors of the structure. Gray arrows are the rung vectors of 
the two types of ladders di,2. (a) 3D structure, (b-d) Projection of the 3D structure onto 
the plane perpendicular to the a, b and c axes. [Taken from Ref. l38]] 
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CHAPTER 3 



Methods 



Due to the exponential growth of the dimension of the Hilbert space with the number 
of sites, the theoretical investigation of quantum many body systems requires highly 
sophisticated techniques. In order to deal with weakly coupled spin- 1/2 ladders, we 
focus here on methods suited for the treatment of quasi one-dimensional systems and 
their mean field extension such as the density matrix renormalization group (DMRG) 
and the Luttinger liquid theory (LL) . 

First, we give an overview of the so called density matrix renormalization group 
(DMRG) or matrix product state (MPS) method. This numerical method introduced 
by S. R. White in the beginning of the 90's is a very powerful technique to treat 
quantum many body physics in particular for one-dimensional systems. The DMRG 
allows one to investigate both static and dynamic properties at zero and finite tem- 
perature. 

Second, we introduce an analytical low energy description for the gapless regime, 
the Luttinger liquid theory (LL). This quantum field theory is the cornerstone of the 
analytical description of many one-dimensional systems. In many situations, the 
bosonization in combination with a numerical determination of its parameters gives 
a quantitative description of the low energy physics. 

As we will see in the following these two methods are complementary and the 
choice depends mainly on the energy or temperature regime we want to focus on. 
Indeed, the DMRG provides a description of the high and intermediate energy prop- 
erties as well as a description of the ground state. However, due to several numerical 
limitations, this method fails to describe the physics at very low energy. In contrast, 
the LL theory focuses on this regime and provides a quantitative description once its 
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parameters are determined for the underlying model (using DMRG for example). 

Finally, we use a mean field approach to treat the weak interladder coupling in the 
case of weakly coupled spin- 1/2 ladders both analytically and numerically. This ap- 
proximation neglects the low quantum fluctuations related to the weak interactions, 
but fully takes into account the fluctuations along the ladders. 



3.1 DMRG 

The DMRG is a numerical method used to determine static and dynamic quantities 
at zero and finite temperature of quasi one-dimensional systems. This method was 
originally introduced by S.R. White [66^ 67 ] to study static properties of one dimen- 
sional systems. Since usually the dimension of the total Hilbert space of a many- 
body quantum system is too large to be treated exactly, the main idea of the DMRG 
algorithm is to describe the important physics using a reduced effective space. This 
corresponds to a variational approach in a space of MPS wave functions. The DMRG 
has been proven very successful in many situations and has been generalized to com- 
pute dynamic properties of quantum systems using different approaches in frequency 
space | [68}j7Q| . Recently the interest in this method even increased after a success- 
ful generalization to time-dependent phenomena and finite temperature situations 
| [39}]4T]|43 - 45J . The real-time calculations further give an alternative route to deter- 
mine dynamic correlation functions of the system | [4Q| which we use in the following. 

In this section, we present an overview of the method providing first a short re- 
view of the basic ideas and focusing on its application to one-dimensional systems 
at zero temperatur^ Furthermore, we discuss the implementation of the time evolu- 
tion and the computation of momentum- frequency correlations. Finally, we extend 
the method to the simulation of finite temperatures. It allows us to compute thermo- 
dynamic quantities as the magnetization or the specific heat. 

More details on the method, its extensions and its successful applications can be 
found in Refs. ^68j-^72J. 



3.1.1 Basic idea of DMRG 

The basic idea of the DMRG technique consists in splitting the Hilbert space in two 



different blocks called system (S) and environment (E) (see Fig. 3.1 ). Each block 



contains a certain number of sites and , respectively. Doing so, a state 

Y,Ca0\w'aW^) (3.1) 

a(3 

of the global Hilbert space (S(^E called superblock) can be decomposed in the basis 

-"I 



formed by the tensor product of the states \w^) and \wq) of the two blocks S and 



^We use in this description a "classic" approach of the DMRG (proposed by S.R. White [67]). 
Overviews over the more flexible MPS formulation can be found in 1 ^68^7 1^^72J . 
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E, respectively. For an exact description of the superblock, it is clear that the size of 
the basis of each block {M^ and M^) would grow exponentially with and . 
Nevertheless, in order to make the computations accessible, it is possible to reduce 
the basis of each block keeping only M < , states to approximate each of 
them. 



S 



OiO 



Figure 3.1: DMRG decomposition of the Hilbert space in two blocks (the system S and 
the environment E). The two circles symbolize the two sites lying at the boundary of 
each blocks for which the description is exact. 

A powerful optimization of the truncated bases is provided by the singular value 
decomposition (SVD) in Ref. |73|. In linear algebra, the SVD of a matrix C G 
Mn,m(C) having the elements Ca(3 ( |3.1| ) is a factorization of the form: 

C = UWV (3.2) 

where W G M^,^(R+) is upper diagonal, and U G M^(C) and V G M^(C) are 
unitary. We denote Wai3 the elements of W with Waa = > and Wa/s = when 
a /3. The elements of U and V are denoted Ua/s and Vais, respectively. Using the 
SVD ( [3.2| ) in order to form a new basis 

of each block, the description of the state {ip) ( |3.1| ) simplifies to the so-called Schmidt 
decomposition |[73|: 

\i^) = J2wa\w^w^) (3.4) 

a 

where Wa are positive numbers with the normalization property = 1 (so 

{ip\ip) = 1). This decomposition is easily computed through the reduced density 
matrices and : 

= TvE\i^){^P\ = J2wl\w^J{w^J 

a 

p^ = Trs\^){^\ = Y,wl\w^){w^\. (3.5) 

a 

The second equality is directly deduced from Eq. ( |3.4| ) and shows that the eigenbases 
\w^) and \ w^) of p^ and p^ correspond to the basis \ w^w^) of the Schmidt decom- 
position ( |3.4| ) respectively. Thus according to Eq. ( |3.5| ), the Schmidt decomposition 
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can be obtained by computing and diagonalizing the reduced density matrices 
and . 

Thereby we optimize the basis of each block for the description of | V^) truncating 
the Schmidt decomposition ( |3.4| ) keeping only the states \ w^) and \ w^) correspond- 
ing to the M largest ' 

M 

W-^Wa\w^W^) = \ll^) with UJi>UJ2>UJ3>"->0. (3.6) 

The quality of this basis for the description of |?/;) can be quantified by the so-called 
truncation error: 

eM= E < = m-m' (3.7) 

which measures the difference between the exact state lip) and its approximation 
\ip) using the optimized truncated basis. For a given M, this error is minimized by 
the truncation procedure ( |3.6| ) and clearly depends on the distribution o^ Wa • In 
the systems treated in this work the DMRG is essentially exact (cm 0) for a 
reasonable choice of M. 



3.1.2 Finite DMRG 

In order to deal with finite dimensional systems, we present the standsLvd finite DMRG 
algorithm that allows one to compute the ground state |?/^o) with open boundary con- 
ditions. 

This method consists in optimizing iteratively the bases \ w^) and \ w^) of dimen- 
sion M for each decomposition of the superblock following the truncation procedure 



described in Sec. 3.1.1 At each step of the method the edge of each block is shifted 



by one lattice site as shown in Fig. 3.2 Next, |?/^o) is computed in the actual two 



block decomposition using a Lanczos method | 74 | a nd the corresponding optimized 



bases \w^) and \ w^) are determined following Sec. 3.1.1 The approximate ground 
state iV^o) in this truncated basis ( [3.6| ) is then kept as the input state in the Lanczos 
method for the next step of the finite DMRG[^ To reach an optimal description of 
IV^o), before converging, this procedure has to be repeated a few times |68 | passing 
through all the sites of the system. The observables are computed during the last 
DMRG sweep when they involve one of the two exactly known sites at the edge of 



the two block decomposition (see Fig. 3.2). Once generated the operators are re 



defined at each step of the DMRG according to the new optimal basis before being 
combined and their expectation value computed. 



^The distribution of the Wa depends on the entanglement of |^} between the two blocks. Although a 
part of this entanglement is lost during the basis truncation \3.6) the latter procedure minimizes this loss 
for the Schmidt decomposition ( |3.4) . 

^\tpo) is considered as good approximation of\ipo) in the next two block decomposition of the system. 
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finite DMRG: 



OiO 



E 



OiO 



OiO 



OiO 



infinite DMRG: 



oo 

]♦ t[ 

oo 

It 



Figure 3.2: Sketch of the finite and infinite DMRG procedures. The two circles sym- 
boHze the two sites lying at the boundary of each block for which the description is 
exact. 



A good starting state for the finite DMRG algorithm is provided by the infinite 
DMRG procedure. This consists in growing the system by introducing additionnal 
sites, two by two, in its center (as shown in Fig. |3.2| ). This recursive procedure 
starts with a small enough system for which the ground state |V^o) can be computed 
exactly. Then the system is decomposed in two symmetrical blocks and their optimal 
M states basis \w^) and \w^) with respect to {ipo) is determined as described in 



Sec. 3.1.1 These bases are used to approximate each block when adding the two 
next new sites at their edge. This choice of basis assumes that the two additional 
sites in the center do not change too much the ground state. This assumption is true 
while the thermodynamic limit is reached, but can be very bad at the beginning of 
the procedure, when the system is small. Nevertheless the infinite DMRG provides 
a usually good starting state for the finite DMRG. 



3.1.3 Time dependent DMRG 

The t-DMRG |39}g2| (time dependent DMRG) method is based on the principle of 
the original DMRG (see Sec. |3.1.1| ). In order to deal with the time simulation, an 
effective reduced Hilbert space is chosen at each time step to describe the physics 
one is interested in. The implementation of this idea can be performed using different 
time-evolution algorithms. Here we use the second order Trotter- Suzuki expansion 



19 



3. METHODS 



for the time-evolution operator of a time-step St |4Q||41| : 



I odd I even I odd 



where i^/ is the local Hamiltonian on the bond linking the sites / and / + 1. This 
decomposition is valid when the total Hamiltonian H decomposes in a sum of such 
terms Hi . 

In order to apply the time-evolution ( |3.8| ) on a given state we use the sweep 



procedure presented in Sec. |3.1.2] for the finite DMRG. During the sweep, each term 
of ( [3.8| ) is applied on the system while it involves the two exactly known sites of the 
two block decomposition (see Fig.|3.2 [). Aft er each one of these steps, the optimized 
basis is updated as described in Sec, p.l.l] Thus the t-DMRG adapts its effective 
description at each time-step. 



In addition to the truncation error discussed in Sec. |3.1.1] present in all variational 
methods, the t-DMRG is also limited by the expansion of the time-evolution opera- 
tor ( |3.8| ). This second uncertainty can be controlled by the choice of the time- step 5t 
(see Ref. |75| for a detailed discussion). 



3.1.3.1 Momentum-energy correlations 

To obtain the spectral functions (Sq^{q,uj) in ( |5.1| )), we first compute the correla- 
tions in space and time 

^^f (^n) = (0|e^^-^^r+L/2,.e-^'^''^f/2,ilO) (3-9) 

with / = -L/2 + 1, -L/2 + 2, . . . , L/2, A: = 1, 2, and tn = nSt (n = 0,1, ... , Nt) 
is the discrete time used. These correlations are calculated by time-evolving the 
ground state |V^o) = |0) and the excited state = S^^^ il^o) using the t-DMRG 
(see Sec. 3.1.3). Afterwards the overlap of S^_^j^^2 fel^i(^)) IV^o(^)) is evaluated 
to obtain the correlation function ( [3.9| ). 

In an infinite system reflection symmetry would be fulfilled. To minimize the 
finite system corrections and to recover the reflection symmetry of the correlations, 
we average them 



(3.10) 



We then compute the symmetric (antisymmetric) correlations (upon leg exchange) 
(see Sec. |5.1[ > 

sriiu) = 2(5,"f (t„) ± sriiu)) (3.11) 

with the rung momentum Qy = 0, tt, respectively. Finally, we perform a numerical 
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Fourier transfornQ 



Nt L/2 



(3.12) 



n=-Nt + l Z=-L/2+l 



for discrete momenta q = 2iTk/L (/c = 0, 1, . . . , L — 1) and frequencies uj. The 
momentum q has the reciprocal units of the interrung spacing a (a = 1 is used if not 
mentioned otherwise). Due to the finite time step 5t, our computed Sq^{q^uj) are 
limited to the frequencies from —7r/Stt0 7r/6t. The finite calculation time tf = NtSt 
induces artificial oscillations of frequency 27r/t/ in Sq^{q^uj). To eliminate these 
artefacts and reduce the effects of the finite system length, we apply a filter to the 
time- space correlations before the numerical Fourrier transform ( |3.12| ), i.e. 



if not stated otherwise. As the effect of this filtering on the momentum-energy cor- 
relations consists to convolve them by a Gaussian function 



it minimizes the numerical artefacts but further reduces the momentum-frequency 
resolution. 

After checking the convergence, typical values we used in the simulations of 
dynamic quantities are system lengths of up to L = 160 sites while keeping a few 
hundred DMRG states M. We limited the final time tf to be smaller than the time 
necessary for the excitations to reach the boundaries (tf ~ L/2u with u the LL 
velocity in Fig. |3.4| ) in order to minimize the boundary effects. The computations 
for the BPCB parameters, Eq. ( |2.15| ), were typically done with a time step of St = 
0.0355 J|| ~^ up to = 71 J|| ~ (but calculating the correlations only every second 
time step). The momentum-frequency limitations are then Suj ^ 0.11 Jy and Sq ^ 
0.1 a~^. Concerning the other couplings and the spin chain calculations, we used a 
time step St = 0.1 Jy"^ up to tf = 100 Jy"^ (also with the correlation evaluations 

^The negative time correlations in the sum for n < are deduced from their value at positive time 

since = S^f (q,tn), with S^^^ (q,tn) = e-^'?^^;^^^ (tn), for translation invariant 

systems, and for correlations such as S*^^ = S^^ with a, (3 = +, — . In order to delete the numer- 
ical artefacts appearing in the zero frequency component of SQ^(q,uj) due to the boundary effects and 
the limitation in the numerical precision, we compute the Fourier transform of S'^^g(tn) only with the 
imaginary part that has no zero frequency component as proposed in Ref. \T6j. 




(3.15) 
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every second time steps) for a momentum-frequency precision 5uj ^ 0.08 J\\ and 
5q ^ 0.1 a-^. 

Different techniques of extrapolation in time (using linear prediction or fitting the 
long time evolution with a guessed asymptotic form cf. Refs. [77. 78|) were recently 
used to improve the frequency resolution of the computed correlations. Neverthe- 
less, as none of them can be applied systematically for our ladder system due to the 
presence of the high energy triplet excitations (which result in a superposition of 
very high frequency oscillations), we decided not to use them. 



3.1.4 Finite temperature DMRG 

The main idea of the finite temperature DMRG | [43] - [45| (T-DMRG) is to represent 
the density matrix of the physical state as a pure state in an artificially enlarged 
Hilbert space. The auxiliary system is constructed by simply doubling the physical 
system. Doing so, we define the totally mixed Bell state on the auxiliary system as 

iV'(o)) = ^nEi^r^F) (3.16) 



where \afaf^) is the state at the bond / of the auxiliary system which has the same 
value \af) = |crp) = |cr^) on the two sites of the bond (the physical (P) and its copy 
(C)). The sum is done on all these A^^- states \ai). Considering the propert}r] 

IVc|^(0))(^(0)| = -fj^ (3.17) 

of the state |V^(0)), it is possible to construct the Boltzmann distribution at finite 
temperature T = 1//3. 

Starting from the infinite temperature limit, we evolve down in imaginary time 
the physical part of \ip{0)) to obtain 

= e-^"/^\m) (3.18) 



using the t-DMRG algorithm presented in Sec. |3. 1.3] with imaginary time. In order 
to avoid an overflow error, this state is renormalized at each step of the imaginary 
time evolution: 

^ MP)) (3 19) 

Hence, according to ( |3.17| ), we get the Boltzmann distribution through 



is the identity operator on the physical system and Trc is the partial trace on the copy system. 
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Therefore, the expectation value of an operator O acting in the physical system with 
respect to the normalized state is directly related to its thermodynamic aver- 

age, i.e. 

{o)p = ^^^f^ = wmm)- (3.21) 

We use this method to compute the average value of the local rung magnetization 
m^{T) and energy per rung E{T) in the center of the system. Additionally we 
extract the specific heat c(T) by 

cW + 6/3/2) ^ {{E)p+sp - {E)p) (3.22) 

where 5P is the imaginary time-step used in the T-DMRG. 

To reach very low temperatures T ^ for the specific heat, we approximate the 
energy by its expansion in T 

n 

E{T)^Eo^^a,T' (3.23) 

i=2 

up to n = 4. The energy at zero temperature Eq is determined by the original 
DMRG in Sec. ( |3.1.1| ). Since E{T) has a minimum at T = the linear term in the 
expansion ( |3.23| ) does not exist. The numbers (z = 2, 3, 4) are obtained fitting the 
expansion on the low T values of the numerically computed E{T). 

After checking the convergence, typical system lengths used for the finite tem- 
perature calculations are L = 80 (L = 100 for the spin chain mapping) keeping a 
few hundred DMRG states M and choosing a temperature step of (5/3 = 0.02 
(^/3 = 0.01 for the spin chain mapping). 

Let us note that recently a new method has been developed to treat finite temper- 
atures |l7T][79]|80l which is very promising to reach even lower temperatures. 



3.2 Luttinger Liquid (LL) 

In quantum systems, the interactions between particles can lead to very different 
physics which depend strongly on the dimensionality of the system. For instance, 
in high dimensions many systems enter into the universality class of Fermi liquids. 
This theory describes systems in which the elementary excitations are quasiparti- 
cles | [8T| . Contrarily in ID systems the effects of the interactions is so strong that 
the excitations are generally collective. The Luttinger liquid theory describes such 
systems in which the collective excitations are free bosonic excitations with linear 
spectrum. In this situation, the physics is described by the Hamiltonian fWf\ 



Hll = ^J dr [uK{drO{r)f + -| {drHr)f 



(3.24) 
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where (j) and 6 are canonically commuting bosonic fields, 

[(^(r), dr'Oir')] = i7rS{r - r'). (3.25) 

Many gapless ID interacting quantum systems belong to the Luttinger Liquid (LL) 
universality class: the dynamics of their low-energy excitations is governed by the 
Hamiltonian ( |3.24| ) and the local operators of the underlying model are written through 
the free boson fields and (j) (the latter procedure is often called bosonization). 

The dimensionless parameter K entering Eq. ( [3.24| ) is customarily called the 
Luttinger parameter, and u is the propagation velocity of the bosonic excitations (ve- 
locity of sound). These parameters are non-universal and depend strongly on the 
underlying model. Once computed all the time and space correlations can be deter- 
mined asymptotically by the field theory corresponding to the Hamiltonian Hll- 

From an experimental point of view, theoretical predictions of the LL theory 
have been observed in a growing number of ID systems such as the organic conduc- 
tors | [82| , quantum wires |83|, carbon nanotubes |84|, edge states of quantum Hall 
effect f85l|, ultracold atoms |37|, antiferromagnetic spin chain |86| or spin ladder 
systems |87|. In these systems characteristic features of the LL theory such as the 



power law behavior of some correlation or spectral functions, discussed in Sees. 3.2.2 



and 3.2.3 for spin-1/2 ladders, have been observed. However, since the details of the 
interactions are rarely known, only a theoretical estimate of the power law exponents 
related to the parameter K was usually possible. 

In this section, we present the LL predictions of the low energy physics of spin- 
1/2 ladders. Computing precisely the parameters of the LL theory for the spin-1/2 
ladder model, we quantitatively describe its low energy properties. Compared to the 
measurements on the compound BPCB in chapter [4j this description provides the 
first quantitative test of the LL theory. 



3.2.1 Bosonization of the spin-1 / 2 chain and ladder 



It has been shown that the gapless regime of the spin-1/2 ladder model ( |2.3| ) is de- 
scribed by the LL theory |[3][2T][2^| . In particular the bosonization of the local spin 
operators performed for both the strong (7 0) and the weak (7 ^ 00) coupling 
limits are smoothly connected |3 1. In the following, we perform the more straightfor- 
ward strong coupling procedure based on the spin chain mapping (Sec.|2.1.l]) starting 



with a reminder of the bosonic description of the mapped spin-1/2 XXZ chain. 

The spin-1/2 XXZ chain, Eq. ( |2.11| ), in the gapless phase is a well-known ex- 
ample of a model belonging to the LL universality class. Its local operators are 
expressed through the boson fields as follows |(6| : 



S^{r) = e 



/2A^(-1)^ + 2V^cos(2(^(r) - 27rmV) 



dr(t){r) 



/2A^(-1)^ cos(2(^(r) - 27rmV). 



(3.26) 
(3.27) 
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Here the continuous coordinate r = la is given in units of the lattice spacing a, 
rh^ = l^M^) /L is the magnetization per site of the spin chain, and Ax, Bx and Az 
are coefficients which depend on the parameters of the model ( [2.11| ). How to calcu- 
late Ax, Bx^ and Az is described in Sec. 3.2.2 For the XXZ spin-1/2 chain, 
a geometrical representation of the two fields and (j) in Eqs. ( |3.26| ) and ( |3.27| ) is 
easily obtained using their classical interpretation. As shown in Fig. |3.3[ they can be 
seen as the two polar angles of the spin fluctuation 5S = S — m^z which gives an 
intuition of the origin of the two terms in the LL Hamiltonian ( |3.24| ). The first term 
(9r^(r))^ measuring the spatial fluctuations of (j){r) is related to the longitudinal (z) 
direction interaction term in ( |2.11| ). In contrast, the transverse (xy) interactions are 
responsible for the second term {drO{r))'^ in ( |3.24| ) related to the spatial fluctuations 
of 0{r). Similarly to the original spin commutation relation ( |2.2| ), the quantum nature 
of the two fields comes from their commutation relation ( |3.25| ). The latter relation 
induces a competition between the ordering in the transverse and the longitudinal 
direction. Due to the strong effects of quantum fluctuations in ID systems, the cor- 
relation functions in a LL decay algebraically with exponents depending on K (see 
Sees. [3:2:21 and |3.2.3.2| ). 




Figure 3.3: Using a classical interpretation, the fields 2(f) and 6 can be viewed as the two 
polar angles of the spin fluctuation SS. 

As discussed in appendix |A| the Hamiltonian ( |2.11| ) is the leading term in the 
strong coupling expansion of the spin-1/2 ladder model ( |2.3| ). Using this strong 
coupling approach, local operators of the latter model are bosonized by combining 
Eqs. ( |2lQ| ), 13:261 ), and ( [3^71 ) 

Sf{r) = {-ly e^^^(^)+^"^ [v^+ \/2^cos(2^(r) - 27rmV) 



(3.28) 

S]{r) = ^ - cos(20(r) - 2^mV). (3.29) 

with j = 1, 2 is the number of the leg. We would like to stress that even for a small 
7 some parameters out of Ax, Bx, and Az show significant numerical differ- 
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ences if calculated within the spin chain ( |2.11| ) compared to the spin ladder ( [2.3| ) (see 
Fig. 3A_). We discuss this issue in Sec. 3.2.2[ 



In the following, we first discuss the numerical determination and the properties 
of the parameters A^, B^, and Az. Furthermore we recall some properties of 
the LL focusing on the finite and zero temperature correlations which are directly 
related to many experimental quantities (see chapters]?] and [5]). 

3.2.2 Luttinger liquid parameter determination 

In this paragraph we detail the determination of the LL parameters u, K and the 
prefactors A^, and A^ (see Eqs. ( [T24| ), ( [^26] ), ( [^27] ), ( [^28] ) and ([T29])) using 
two main properties of the LL ground state i.e. the algebraic decay of the correlation 
function^ ( |3.31| ) and ( |3.32| ) as well as the susceptibility^ ( |3.3Q| ). These parameters 
are necessary for a quantitative use of the LL theory. The parameters A^, B^ 
and Az and their dependence on the magnetic field have been previously determined 
in Refs. \Tlj2^ 88] for different values of the couplings than those considered here. 
We obtain these parameters in two steps |l3T||89l : 



(i) We determine the ratio u/K from its relation to the static LL susceptibility^ jSj 
|90) 

We numerically compute the static susceptibility using DMRG and infer the 
ratio u/K with a negligible error. 

(ii) The parameter K and the prefactors Ax, B^ and Az are extracted by fitting 
numerical results for the static correlation functions obtained by DMRG with 
their analytical LL expression^ [28J 

{bibi,)-Ax- — --^-Bx[-l) - — (3.31) 

These correlations computed for infinite systems decay algebraically with the 
distance |/ — V\ with a K dependent exponent. In practice, we use the more 
sophisticated expressions Eqs. ( |B.2| ), ( |B.3| ) and ( |B.4| ) (for {S[)) shown in 
appendix |B] and derived in Ref. |28|. These take into account the boundary 
effects which are present in the finite DMRG computations but neglected in 
Eqs. ( [33T] ) and ( [332] ). 



^ These relations are given for the spin chain \2.1l) si nce from these the relations for the spin ladders 
can be easily inferred using the spin chain mapping \B.l) . 
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We first fit the transverse correlation (xx -correlation (SfSf,)) to extract the 
parameters K, Ax, and Bx. Then we use the previously extracted value for K 
to fit the longitudinal correlation (zz-correlation (SiSf,)) and the magnetiza- 
tion, {S[), which allow us to determine Az. The values determined by both 
fits are very close and in Fig.|3.4|the average value of both is shown. 



All the results presented in Fig. 3.4 were obtained for L = 200 and several 
hundred DMRG states after an average on the four sets of used data points in the fit 
10 <l,r < 170, 30 <l,r < 170, lO < l,r < 190, 30 <l,r < 190. The error 
bars correspond to the maximum discrepancy of these four fits from the average. We 
further checked that different system lengths lead to similar results. 




0.2 0.4 0.6 0.8 1 



Figure 3.4: LL parameters u, K and the prefactors of the spin operators Ax, Bx, Az 
versus the magnetization per rung, , computed for a spin ladder with the BPCB cou- 
plings \2.\5\ (red crosses) and for the spin chain mapping (blue stars). The strong cou- 
pling expansion of u and K up to second order in 7 (discussed in appendix |A.4| ) is 
plotted in black dashed lines. [Taken from Ref. |38 |] 
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The LL parameters of the ladder system ( |2.3| ) for the BPCB coupHngs ( |2.15| ) are 
presented in Fig. 3.4 as a function of the magnetization per rung. Additionally we 
show the parameters of the spin chain mapping (computed for the spin chain Hamil- 
tonian ( |2.11| )) for comparison. When the ladder is just getting magnetized, or when 
the ladder is almost fully polarized, K ^ 1 (free fermion limit) and u ^ (be- 
cause of the low density of triplons in the first case, and low density of singlets in 
the second case). Between these two limits K < 1 due to the triplet-triplet repulsion 
(see Eq. ( |A.13| ) in the strong coupling expansion). For the spin chain mapping, the 
reflection symmetry around = 0.5 arises from the symmetry under tt rotation 
around the x or y axis of the spin chain. This symmetry has no reason to be present 
in the original ladder model, and is an artefact of the strong coupling limit, when 
truncated to the lowest order term as shown in appendix |A| The values for the spin 
ladder with the compound BPCB parameters can deviate strongly from this sym- 
metry. The velocity u and the prefactor remain very close to the values for the 
spin chain mapping. In contrast, the pref actors Az, Ax and the exponent K deviate 
considerably and A^ and K become strongly asymmetric. The origin of the asym- 
metry lies in the contribution of the higher triplet states pi, and can be understood 
using a strong coupling expansion of the Hamiltonian ( |273| ) up to second order in 
7 (see appendix |A.4| ). This asymmetry has consequences for many experimentally 
relevant quantities and it was found to cause for example strong asymmetries in the 
3D order parameter, its transition temperature and the NMR relaxation rate as will 
be discussed in chapter]?] (see Fig. |4.5[ Fig. |4.6| and Fig.|4.7). 



3.2.3 Dynamical correlations 

In this paragraph, we focus on the retarded correlations defined in the time-space as 



(3.33) 



with 6(t) the Heaviside function, a, /3 = ±, z such as = S^^ and 6'f (r, t) = 
e^^^Sf (r)e~*^^ is the time evolution of the operator Sf (r) . Their Fourier transform 



is computed as {q^^) 



J dr dt 



J{u;t-qr) af3 



x'iji'^^'^)- These correlations are 



necessary for the mean field determination of the transition temperature Tc to the 3D- 
ordered phase (see Sec . 4.4. They are also directly relat ed to the NMR relaxation 
rate T^^ (see Sec. 4.3 ) and the INS cross-section (see Sec. 5.6) through the spectral 
functions 



s: 



cx(3 



/oo 
dr dt e^("*-«'^)(5f(r,t)5f (0,0)). 
-OO 



(3.34) 
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As discussed in Sec . |3.1.3| these spectral functions are also accessible numerically 
at zero temperature ( |3.12| ) and have the following properties 

(3.35) 

Hence, the spectral function S^j^ diverges in the low energy limit 1 ^uj unless 

the correlation xtf vanishes in this limit. At zero temperature, S^j^ vanish for all 

negative frequencies. As we will see in Sec. 5.1 measure the excitations of 
the system, their vanishing thus physically means tnat at zero temperature no excited 
state has an energy lower than the ground state. 



3.2.3.1 Finite temperature LL correlations 

Using the bosonization formalism ( |3.28| ) and ( |3.29| ), and taking into account only the 
most relevant terms, we can compute the Fourrier transform of the correlations ( [3.33| ) 
as described in Ref. |[6| for the LL Hamiltonian ( |3.24| ): 



Xi]^{q.^) =(-l)^+^A, A/4k(^-^,^,/3) (3.36) 

1 

4^' 



1 



(3.37) 



with 



r. / ^^ sin (TTiy) / 27r \ ^ / 0((^ — uq) v , 

/ A ^^^^ /a ao^ 

[yj + i\)^Y — w^q"^ 

where B{x^y) = ^-^{I^^y) • correlations ( |3.36| ) and ( |3.37| ) are linear combina- 
tions of the functions and which depend only on the two LL parameters u and 
K. The weight of the /j, component is related to the pref actors and Az for Xip" 
and xf/, respectively. In contrast the component gj^ in ( [3.37 ) has a constant prefactor 



and is invariant in temperature. As we will see in chapter 4j these correlations are 
necessary to compute the critical temperature Tc of the 3D transition through a mean 
field treatment of the interladder coupling J^ 
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Thus using ( [3.36| ), ( [3.37| ) and ( [3.35| ), the local correlations (momentum average 
of the spectral functions ( |3.34| )) 

S"P{lo) = Sf{r = 0,oj) = ^ J dq S^f{q, co) (3.40) 



become in the low energy limit 1 Puj: 



^ ' u \j3u) \4.K' 2K 



A^cos{'kK) (2t:V^ ^ K 



(3.41) 



In this low energy limit, the local correlations ( [3.41| ) are directly related to the NMR 
relaxation rate T^^ , Eq. ( |4.2| ). 



3.2.3.2 Zero temperature correlations in the LL 



At zero temperature the correlation functions in the LL have been computed in 
Ref. ||2T 22 9TJ. In the following, we give directly the expression of the symmetric 



(+) and antisynmietric (— ) spectral functions S^^ = 2{Si{ ± with rung mo 
mentum = 0, tt, respectivel}Q These are the relevant quantities for a comparison 
with INS measurements (see Sec. |5.6| ). They are derived analogously to the finite 
temperature correlation^ ( |3.36| ) and ( |3.37| ) in the limit /3 — > 00 using Eq. p.35| ): 



S^'iq, uj) = (2iTm'f5{q)5(uj) + ^Q(uj) \5(uj - uq) + 5(uj + uq)] 

u 



uT{KY 



Q{uj -u\q-2iTm^\) 



uj^ — u'^{q — 27rm^)^ 



l-K 



(3.43) 



^Note that the edge exponents in the incommensurate branches of the correlations \3A4] are inverted 
compared to their expression in Ref. p2j and pictured in Fig. |3.5| b-c. 

^Note that the two last terms of the zero temperature spectral functions originate from those 

which mix the two fields and cj) in the bosonic description. These were neglected in the bosonic deriva- 
tion of the finite temperature correlations j3.36) . 
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2u 



9(cj — u\q — 7r|) 



e{uj - u\q - 7r{l - 2m^)\) 



1-1/4K 



^00 — u[q — 7r(l — 2m^)] 



2u 



- u[q — 7r(l — 2m^ 



1-77+ 



(3.44) 



with r]± = l/AK ±1 -\- K. The spectral function + is obtained replacing 
—m^ in the expression Eq. ( [3.44| ). 



> (a) 




7t{1 - 2m^) TT 7r(l + 2m^) 




7r(l - 2m^) TT 7r(l + 2m^) 



Figure 3.5: Map of the low energy spectral functions of the LL model where the white 
areas represent the continuum of excitations. In the striped areas no excitations are 
possible, (a) Sq^ {q, uj): the dash-dotted lines (blue) are the excitation peaks close to ^ = 
0, 27T and the dashed lines (red) are the continuum lower boundary with edge exponent 
1 — K close to ^ = 27rm^, 27r(l — m^). (b) {q, cj), (c) S~~^{q, cj): the continuum 
lower boundary close to ^ = tt, 7r(l di 2m^) is represented by solid lines (black) (edge 
exponent 1 — 1/4K), dashed lines (red) (edge exponent 1 — rj- — 2 — 1/4K — K) and 
dash-dotted lines (blue) (edge exponent 1 — ?7+ = — 1/ 4K — K). [Taken from Ref. 1 38 1] 



The expressions Eq. ( |3.43| ) and Eq. p.44| ) exhibit the typical behavior of the 
frequency-momentum LL correlations: a continuum of low energy excitations exists 
with a linear dispersion with a slope given by the Luttinger velocity ±.u. The spectral 
weight at the lower boundary of the continuum displays an algebraic singularity with 
the exponents related to the Luttinger parameter K. A summary of this behavior is 



sketched in Fig. 3.5 For the considered system the longitudinal correlation S^^ is 



predicted to diverge with the exponent 1 — K at its lower edge. As shown in Fig. 3.6 b 
the exponent of this divergence is very weak < 0.2 for the parameters of BCPB. The 
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Figure 3.6: Different exponents that appear in the LL correlation functions, Eqs. ( |3.43| ) 
and ( |3.44| ), versus the magnetization . The soHd (dashed) Unes are determined from 
the ladder (spin chain mapping) exponent K in Fig. |3.4| The exponent 1 — K of the Sq^ 
correlations is shown in (b), and the exponent 1 — l/AK of the S^^ correlations at the 
q = 7T branch in (a) (lower red curves). The exponents l—rj- —2 — 1/4K — K (upper 
black curves) in (a) and 1 — ?7+ = — 1 / — K in (c) correspond to both sides of the 
incommensurate branches of the (see Fig. 3.5 1. [Data taken from Ref. L38J] 



transverse correlations S^^ exhibit a distinct behavior depending on the considered 
soft mode. Close to g' = tt the weight diverges with an exponent given by 1 — 1/AK. 
This divergence is strong for the considered parameters (1 — 1/4:K ^ 3/4 ^ in 
Fig. 3.6 a). In contrast at the soft mode q = 7r(l — 2m^), 7r(l + 2m^) a divergence 
(cusp) is predicted at the lower edge with the exponent 2 — 1 / 4:K — K ^ 3/4 in 
Fig.[l6]a (-1/4K -K^ -5/4 in Fig.[l6]c). 



3.3 Mean field approximation 

Up to now, we have presented methods adapted to deal with one dimensional sys- 
tems. In real compounds, an interladder coupling is often present. As discussed in 



Sec. 2.1.2 in the incommensurate regime this interladder coupling J' (cf. Eq. \2.\) ) 



can lead to a new three dimensional order (3D-ordered phase in Fig. 2.2 b) at temper- 



atures of the order of the coupling J' . In the case of BPCB the interladder coupling 



is much smaller than the coupling inside the ladders, i.e. J' <C J±, J\\ (Sec. 2.2). 
Therefore, unless one is extremely close to hd or hc2 one can treat the interladder 
coupling within a mean field approximation. Let us emphasize that this approach 
incorporates all the fluctuations inside a ladder. However, it overestimates the effect 
of J' by neglecting quantum fluctuations between different ladders. Such effects can 
partly be taken into account by a suitable change of the interladder coupling f32| to 
an effective value that will be discussed in Sec. 14.41 Close to the critical fields the 
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interladder coupling J' becomes larger than the effective energy of the one dimen- 
sional system. This forces one to consider a three dimensional approach from the 
start and brings the physics of the system in the universality class of Bose-Einstein 
condensation |[3][9|. In the following we consider that we are far enough (i.e. by an 
energy of the order of J') away from the critical points so that we can use the mean 
field approximation. 

The mean field approximation of the interladder interactions in the 3D Hamilto- 
nian i^3D (Eq. \2A\ ) reads 

• ^l',k',^i' — ^l,k,fi • (S^/,/c',^') + {^l,k,fi) • — {^l,k,fi){^U,k',fi') (3.45) 

and the ladders decouple. Since the single ladder correlation functions along the 
magnetic field direction (z axis) decay faster than the staggered part of the ones in 
the perpendicular xy plane (see Eqs. ( |3.31| ) and ( |3.32| ) for the LL exponent K of 



the ladder shown in Fig. 3.4), the three dimensional order will first occur in this 
plane. Thus the dominant order parameter is the q = n staggered magnetization 
perpendicular to the applied magnetic field. Focusing on one of the ladders /i of the 
system, we thus introduce the order parameters 

(5,^,) = -(-l)'+'=mS and (5^,) = ^ - (-l)'+'=m^ (3.46) 

assuming the staggered xy ordering to be along the x axis, will be very small 
and therefore neglected. Hence H^j) ( |2.1| ) becomes 

Hu, = J||F|| + J^H^ + Slk + E(-l)'^'^^fc- (3.47) 

l,k l,k 

Here we suppose that the coupling is dominated by Uc neighboring ladders which are 
antiferromagnetically ordered (along the x axis) with respect to each other, where Uc 



is the rung connectivity {ric = 4 for the case of BPCB, cf. Fig. 2.6 ). This mean field 
Hamiltonian corresponds to a single ladder in a site dependent magnetic field with a 
uniform component in the z direction and a staggered component in the x direction. 
The ground state wave function of the Hamiltonian must be determined fulfilling the 
self-consistency condition for and using numerical or analytical methods. 



3.3.1 Numerical mean field 

The order parameters and can be computed numerically by treating the mean 
field Hamiltonian Hmf self consistently with DMRG. These parameters are evalu- 
ated recursively in the center of the ladder (to minimize the boundary effects) starting 
with = and = 0.5. An accuracy of < 10~^ on these quantities is quickly 
reached after a few recursive iterations (typically ~ 5) of the DMRG keeping few 
hundred DMRG states and treating a system of length L = 150. We verified by 
keeping as well the alternating part of the z order parameter that this term is 
negligible (< 10"^). 
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3.3.2 Analytical mean field 



Using the low energy LL description of our ladder system (see Sec. 3.2), it is pos- 
sible to treat the mean field Hamiltonian Hmf within the bosonization technique. 
Introducing the LL operators ( |3.28| ) and ( |3.29| ) in Hmf ( |3.47| ) and keeping only the 
most relevant terms leads to the Hamiltonian [92(|93| 



SG 



1 

2^ 



dr 



uK {drO{r)Y 



dr cos{0{r)) 
(3.48) 

where we neglected the mean field renormalization of in ( |3.47| ). This Hamiltonian 
differs from the standard LL Hamiltonian i^LL ( |3.24| ) by a cosine term corresponding 
to the X staggered magnetic field in ( [3.47| ). It is known as the sine-Gordon Hamil- 
tonian |6, 9£,25|. In the range of the typical K values for BPCB (see Fig. |3.4| ) 
the cosine term in ( |3.48| ) is relevant |[6j and orders the field 0{r). As pictured in 
Fig. [33] this ordering is responsible for the staggered transverse magnetization m^. 
The expectation values of the fields can be derived exactly using integrability [96J . 
In particular can be determined self-consistently as shown in Sec. 4.4.2 
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CHAPTER 4 



Static properties and NMR relaxation rate 



In chapter [2j we have seen that the physics of weakly coupled spin- 1/2 ladders is 
particularly rich. In the following, we explore the diversity of their phase diagram, 
pictured in Fig. |2.2[ by computing several physical quantities such as the magneti- 
zation, the rung state density and the specific heat. In particular, we test the LL low 
energy prediction and evaluate the related crossover to the quantum critical regime. 
Furthermore we discuss the effect of the 3D interladder coupling computing the stag- 
gered magnetization in the 3D-ordered phase and its critical temperature. We finally 
discuss the NMR relaxation rate in the LL gapless regime related to the low energy 
dynamics. All of these physical quantities are computed for the BPCB parameters 
(see Sec. |2.2| ). Hence they can be directly compared to the experiments on BPCB 
discussed in detail at the end of this chapter. 



4.1 Critical fields 

The zero temperature magnetization contains extremely useful information. Its be- 
havior directly gives the critical values of the magnetic fields hd and hc2 at which 



the system enters and leaves the gapless regime, respectively (Fig. 2.2 b). In Fig. 4.1 
the dependence of the magnetization on the applied magnetic field is shown for a 
single ladder and for weakly coupled ladders. At low magnetic field, < hd, the 
system is in the gapped spin liquid regime with zero magnetization, and spin singlets 
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on the rungs dominate the behavior of the systen£] see Fig. 4.2 At = h d, th e 
Zeemann interaction closes the spin gap to the rung triplet band (Fig. 2.2). 
Above > hd the triplet band starts to be populated leading to an in- 
crease of the magnetization with h^. The lower critical field in a 13th order ex- 
pansion |T2| in 7 is /ici ^ 6.73 T for the BPCB parameters. At the same time the 
singlet and the high energy triplets occupation decreases as shown in Fig. |4.2[ For 
> hc2 = J± + 2J|| ^ 13.79 T (for the compound BPCB), the band is 
completely filled and the other bands are depopulated. The system becomes fully 
polarized (m^ = 1) and gapped. The two critical fields, hd and hc2, are closely 
related to the two ladder exchange couplings, J± and Jy. As they are experimen- 
tally easily accessible, assuming that a ladder Hamiltonian is an accurate description 
of the experimental system, these critical fields can be used to determine the ladder 
couplings |31|. 

Such a general behavior of the magnetization is seen for both the single lad- 



der and the weakly coupled ladders in Fig. |4.1| In particular, the effect of a small 
coupling J' between the ladders is completely negligible in the central part of the 
curve. Only in the vicinity of the critical fields, the single ladder and the coupled lad- 
ders show a distinct behavior. The single ladder behaves like an empty (filled) one- 
dimensional system of non-interacting fermions which leads to a square-root behav- 
ior (X {h^ — hdY^'^ close to the lower critical field and 1 — ex (hc2 — h^Y^'^ 
close to the upper critical field. In contrast, in the system of weakly coupled lad- 
ders, a 3D-ordered phase appears at low enough temperatures in the gapless regime 



(see Sees. 2.1.2 and 3.3). The magnetization dependence close to the critical fields 
becomes linear, oc — h^P, and 1 — oc h^^2 ~ respectively l[3j[54|. 
In comparison with the single ladder, the critical fields h^^^ and h^^2 shifted by 
a value of the order of J' in comparison with hd and hc2 • This behavior is in the 
universality class of the Bose-Einstein condensation |3,9|. Appearing very close to 
the critical fields these 3D effects are at the limit of validity of the mean field ap- 
proximation. Nevertheless they are qualitatively reproduced by this approximation 
as shown in the insets of Fig. |4.1| 

For comparison, the magnetization of a single ladder in the spin chain mapping 



is also plotted in Fig. |4.1[ This approximation reproduces well the general behavior 
of the ladder magnetization discussed above. However, note that for the exchange 
coupling constants considered here the lower critical field in this approximation is 
different from the ladder one. The lower critical field is h^^^^ = J± — J\\ ~ 
6.34 T < hd- The upper critical field h^2^'^ = Jx + 2 Jy = hc2 is the same as 
for the ladder. If we rescale h^^^ and h^2^^ to match the critical fields hd and hc2 

(h^ "^xxzL^j^xx^^""^^ + hd)^ the magnetization curve gets very close to the 

one calculated for a ladder. However, in contrast to the magnetization curve for the 
ladder, the corresponding curve in the spin chain mapping is symmetric with respect 



^The perturbative expression of the ground state in the spin Hquid regime and the corresponding singlet 
and triplet densities are given in appendix [A| 
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Figure 4.1: Dependence of the magnetization per rung on the magnetic field at 
zero temperature for the single ladder computed by DMRG with the BPCB couplings 
(Sec. |2.2| ) (dashed red line), the spin chain mapping (dotted blue line) rescaled to fit with 
the single ladder critical fields (dash-dotted blue line), and for the weakly coupled lad- 
ders treated by the mean field approximation (solid black line). The insets emphasize 
the different behavior of the magnetization curves for the single (dashed red line) and 
weakly coupled (solid black line) ladders close to the critical fields which are indistigu- 
ishable in the main part of the figure. The dotted lines in the insets correspond to the 
linear and square root like critical behavior. See also Fig |4.8| for a comparison with the 
NMR measurements on BPCB. [Taken from Ref. (SSj] 

to its center at h-^^ = ^ = J±-\-J\\/2 due to the absence of the high energy 
triplets. 

4.2 The Luttinger liquid regime and its crossover to 
the critical regime 

The thermodynamics of the spin-1 / 2 ladders has been studied in the past for different 
regimes and coupling ratio 7 [13|23f[25|33| . We here summarize the main interesting 
features of the magnetization and the specific heat focusing on the crossover between 
the LL regime and the quantum critical region using the BPCB parameters (Sec. |2.2| ). 
As the interladder exchange coupling J' is supposed very small compared to the 
ladder exchange couplings Jy and J±, it is reasonable to neglect J' in the regime far 
from the 3D phase. Therefore we focus on a single ladder in the following. 

4.2.1 Finite temperature magnetization 

We start the description of the temperature dependence of the magnetization, m^, 
in the two gapped regimes: the spin liquid phase and the fully polarized phase. For 
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Figure 4.2: Rung state density versus the applied magnetic field at zero temperature 
for the single ladder computed by DMRG with the BPCB couplings. The dash-dotted 
(black) lines correspond to the singlet density (ps). The triplet densities are represented 
by the solid (red) lines for (yo+), the dashed (blue) lines for (po) and the dotted (green) 
lines for (p-). [Taken from Ref. (38]] 



small magnetic fields < hd, the magnetization vanis hes e xponentially, m 



4.3 



a, this decay slowly 



y/7^g-(^ci-/i^)/T^ at low temperature. As shown in Fig. 
disappears while the gap hd — is closed (h^ hd). After a maximum at 
intermediate temperatures decreases to zero for large temperatures due to strong 
thermal fluctuations. Similar features appear for large magnetic fields > hc2- As 
shown in Fig. |4.3[ c, the magnetization increases exponentially up to = 1 at low 
temperature, 1 — ex V^e"'^^^"^^^)/^^ and decreases monotonously in the limit 
of infinite temperature. As in the spin liquid phase, the low temperature exponential 
behavior becomes more pronounced while the gap — hc2 increases. 

In the gapless regime, the magnetization at low temperature has a non-trivial be- 



havior that strongly depends on the applied magnetic field. As shown in Fig. 4.3 b, 
in this regime (hd < < hc2) new extrema appear in the magnetization at low 
temperature. This behavior can be understood close to the critical fields where the 
ladder can be described by a one-dimensional fermion model with negligible inter- 
action between fermions. Indeed, in this simplified picture |6| and in more refined 
calculations l|23j[24][^ the magnetization has an extremum where the temperature 
reaches the chemical potential, i.e., at the temperature at which the energy of ex- 
citations starts to feel the curvature of the energy band. The type of the low tem- 
perature extrema depends on the magnetic field derivative of the LL velocity |[26| 
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(du/dh^). Thus a maximum (minimum) is expected if du/dh^ < (du/dh^ > 0). 
This specific behavior is illustrated in Fig. |4.3[ b with the curve for = 11 T 
(hm = ^^^^4^ < < hc2) with {du/dh^)\h^^iiT < (see Fig. |I4). The 
low temperature maximum moves to higher temperature for < hm and crosses 
over to the already discussed maximum for /i^ < hd (see Fig. 43 a). Symmetrically 
with respect to /i^, a low temperature minimum appears in the curve for = 9 T 
(hci < < hm) with {du/ dh^)\h^ ^9T > (see Fig. 3.4). This minimum slowly 
disappears when h^ hm for which {du/dh^)\h^ ^ (the curve for h^ = lOT is 
close to that). 

The location of the lowest extremum is a reasonable criterion to characterize the 
crossover temperature between the LL and the quantum critical regime |[23][24][26|, 
since the extremum occurs at temperatures of the order of the chemical potential. A 
plot of this crossover temperature versus the magnetic field is presented in Fig. |4.3| e. 
Following this criterium, the crossover has a continuous shape far from hm • Never- 
theless, close to hm we have du/dh^ ^ and the low energy extremum disapears. 
The criterium is thus not well defined and presents a discontinuity at hm which is 
obviously an artefact. In the vicinity of hm, we thus use another crossover criterium 
based on the specific heat (see Sec. |4.2.2] ) that seems to give a more accurate descrip- 
tion. 

The temperature dependence of the magnetization of the spin chain mapping, 
Fig. 4.3 d, exhibits a single low temperature maximum if h^^ < h^ < h^f^ (min- 
imum if h^f^ < h^ < h^^). The appearance of a single extremum and its conver- 
gence to 0.5 when T ^ oo is due to the exact symmetry with respect to the 
magnetic field h^^. This approximation reproduces the main low energy features 
of the ladder but fails to describe the high energy behavior which strongly depends 
on the high energy triplets. 



4.2.2 Specific heat 



Similarly to the magnetization discussed in Sec. 4.2.1 the specific heat of spin-1/2 
ladders shows in the spin liquid and fully polarized phases the typical behavior of 
gapped regimes. At low temperature the specific heat grows exponentially: c ex 

2--3/2g-(^ei-^")/T ^^^^ ^ rj.-2,/2^-{h^-h,2)/T for both phases respectively. After 

reaching a maximum when the gapped excitations start to be thermally populated. 



in the quantum critical regime (see Fig. |2.2| b), it slowly decreases to zero at high 
temperature due to the strong temperature fluctuations. These specific features are 
shown in Figs. 4.4 a and 4.4 c where c(T) is plotted for various applied magnetic 



fields in both gapped regimes. 



As presented in Fig. 4.4 b, the behavior of the specific heat becomes more subtle 
in the gapless regime in which the contribution due to the gapless spinon excitations 
appear at low temperature. This results in a peak around T ~ 1.5 K. This peak is 
most pronounced for the magnetic field values lying mid value between the two crit- 
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Figure 4.3: Temperature dependence of the magnetization per rung, m^(T), for the 
ladder with the BPCB couphngs \2.15\ (a) in the spin Uquid regime for /i^ = 2 T (soUd 
blue lines), — AT (dash-dotted green lines) and /i^ = 6 T (dashed red lines), (b) 
(and (d) for the spin chain mapping) in the gapless LL regime for /i^ = 9 T (solid blue 
lines), /i^ = 10 T (dash-dotted green lines), — 11 T (dashed red lines) and (c) in 
the fully polarized regime for /i^ = 15 T (solid blue lines), — 17 T (dash-dotted 
green lines) and = 19 T (dashed red lines). The results were obtained using T- 
DMRG. The stars at T = K are the ground state magnetization per rung determined 
by zero temperature DMRG. The triangles (squares) mark the low (high) temperature 
extrema. (e) Crossover temperature Tll of the LL to the quantum critical regime versus 
the applied magnetic field (blue circles for the extremum in (T) | criterium and red 
crosses for the maximum in c(T)|^^ criterium). See also Fig. 4.11 for a comparison 



with magnetocaloric effect measurements on BPCB. [Data taken from Ref. p8|] 



ical fields. At higher temperatures the contribution of the gapped triplet excitations 
leads to a second peak which exists also in the gapped regimes as discussed above 
and shown in Figs. |4.4| a and |4.4| c. Its position depends on the magnetic field (see 
Ref. |[23 1 for a detailed discussion). To separate out the contribution from the low 
lying spinon excitations, we compare the specific heat of the ladder to the results 
obtained by the spin chain mapping in which we just keep the lowest two modes 
of the ladder (see Sec. |2.1.l] and appendix |A]). The resulting effective chain model 
is solved using Bethe ansatz |38| and T-DMRG methods. The agreement between 
these methods is excellent and the corresponding curves in Fig. |4.4| b can hardly be 
distinguished. However, a clear difference with the full spin ladder result is revealed. 
While at low temperatures the curves are very close, the first peak in the spin chain 
mapping already lacks some weight, which stems from higher modes of the ladder. 
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In the inset of Fig. |4.4[ b, the low temperature region is analyzed in more detail. 
At very low temperatures the spinon modes of the ladder can be described by the LL 



theory (see Sec. 3.2) which predicts a linear rise with temperature |6^p7j inversely 



proportional to the spinon velocity u (shown in Fig. |3.4| versus the applied magnetic 
field for the BPCB couplings) 



cll(T) 



Ttt 
3u 



(4.1) 



In the inset of Fig. [4!4| b we compare the results of the LL, the Bethe ansatz p8| and 
the DMRG results for the effective spin chain and the numerical DMRG results tak- 
ing the full ladder into account. The numerical results for the adaptive T-DMRG at 
finite temperature are extrapolated to zero temperature by connecting algebraically 
to zero temperature DMRG results (see Sec. |3.1.4| ). A very good agreement be- 
tween ( |4.1| ) and numerics is found for low temperatures. However, at higher temper- 
atures, the slope of the T ^ LL description slightly changes with respect to the 
curves calculated with other methods. This change of slope reflects the fact that the 
curvature of the energy dispersion must be taken into account when computing the 
finite temperature specific heat, and this even when the temperature is quite small 
compared to the effective energy bandwidth of the system. The effective spin chain 
and the numerical results for the ladder agree for higher temperatures (depending on 
the magnetic field), before the higher modes of the ladder cause deviations. 

the location of the low temperature peak 



As for the magnetization (Sec. 4.2.1 



can be interpreted as the crossover of the LL to the quantum critical regime. Indeed, 
in a free fermion description which is accurate close to the critical fields, this peak 
appears at the temperature for which the excitations stem from the bottom of the en- 



ergy band. The corresponding temperature crossover is compared in Fig. |4.3[ e to the 
crossover temperature extracted from the first magnetization extremum (Sec. |4.2.r ). 
The two crossover criteria are complementary due to their domain of validity. The 
first specific heat maximum is well pronounced only in the center of the gapless 
phase. In contrast in this regime the disappearance of low energy extrema renders 
the magnetization criterium very imprecise (cf. Sec. |4. 2.1] ). In Fig. 4.1 1 both criteria 
have been applied on the magnetocaloric effect and specific heat measurements on 
BPCB p3| . These experimentally extracted crossovers are in perfect agreement with 
the ones computed by T-DMRG. 

More generally comparisons of the computed specific heat with actual experi- 
mental data p3| for BPCB are excellent (see Fig. |4.1Q| in the experimental Sec. [43). 
For these comparisons the theoretical data are computed with g = 2.06 related to the 
experimental orientation of the sample with respect to the magnetic field |33 | (see 



Sec. 2.2) and rescaled by a factor 0.98 in agreement with the global experimental 
uncertainties 



^An additional scaling factor 7.47 mJ/gK has to be applied on the theoretical specific heat (per rung) 
to convert to the experimental units. 
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Figure 4.4: Specific heat per rung c versus temperature T computed with T-DMRG for 
the BPCB coupUngs \2A5\ (a) in the spin Uquid regime for /i^ = 2 T (sohd blue fine), 
/i^ = 4 T (dash-dotted green fine), and /i^ = 6 T (dashed red fine), (b) in the gapless 
LL regime for /i^ = 9 T (sohd blue line), = 10 T (dash-dotted green line), and = 
11 T (dashed red line) and (c) in the fully polarized regime for /i^ = 15 T (solid blue 
line), /i^ = 17 T (dash-dotted green line), and /i^ = 19 T (dashed red line). In (b) the 
spin chain mapping at /i^ = 10 T solved by T-DMRG (Bethe ansatz from Ref. 1 38 1) are 
plotted in dashed (dotted) black line. Note that the two lines are hardly distinguishable. 
The triangles (squares) mark the low (high) temperature maxima of the specific heat 
versus temperature. The vertical dashed line marks the temperature T = 0.4 K below 
which the DMRG results are extrapolated (see Sec. |3.1.4| ). The inset in (b) shows the 
low temperature dependence of the specific heat per rung for — 10 T. The T-DMRG 
calculations are in red thick lines for the ladder with the BPCB couplings ( |2.15| l (black 
thick lines for the spin chain mapping). The two curves can hardly be distinguished. 
Their low temperature polynomial extrapolation is plotted in thin lines below T = 0.4 K 
(represented by a vertical dashed line). The linear low temperature behavior of the LL is 
represented by dashed lines (red for the ladder, black for the spin chain mapping). The 
dashed yellow lines correspond to the Bethe ansatz 1 38 1 computation for the spin chain 
mapping. See also Fig. |4.10| for a comparison with measurements on BPCB. [Data taken 
from Ref. |j38J] 

4.3 Spin-lattice relaxation rate 

The NMR spin-lattice relaxation in quantum spin systems is mainly due to the pure 
magnetic coupling 
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between electronic and nuclear spins S = {S^ , S'^ ^ S^) and I = re- 
spectively. Aaf3 with a, /3 = X, z is the hyperfine tensor related to the dipolar 
interaction between the spins S and I. 7n is the nuclear gyromagnetic ratio of the 
measured nuclear spin I (7^ = 19.3 MHz/T for the measurements on BPCB done 



on the N(l) atoms in Fig. 2.5 (see Ref. |31 ])). Due to this form of the coupling the 



NMR spin-lattice relaxation rate is directly related through the Redfield equa- 
tions pSJ to the local correlations S^^{ujo) defined in Eq. ( |3.4Q| ) 

Tf 1 = jlAlS+-{wo) + jlAp^^icoo). (4.2) 

ujo = h^ln is the Larmor frequency. A\\ and A± are the longitudinal and transverse 
components of the hyperfine tensor which have the same order of magnitude than 
the components Aa(3 . 
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Figure 4.5: Magnetic field dependence of the NMR relaxation rate, T^^{h^), at 
T — 250 mK. The solid red line is the bosonization determination using the ladder 
LL parameters for the BPCB couplings shown in Fig. |3.4| (the dashed blue line uses the 
LL parameters of the spin chain mapping). The black circles are the measured NMR 



relaxation rate on BPCB done on the N(l) atoms in Fig. 2.5 from Ref | ,31J . [Taken 
from Ref. (3^] 

Assuming Jy ^ the relaxation rate T^^ can be computed in the gapless 
regime using the LL low energy description (Sec. |3.2| ) of the electronic spin dy- 
namics. We introduce only the most relevant local LL correlation ( |3.41| ) (in the limit 
T > cjo) into Eq. ( |43| ) (i.e. the transverse component (cj ^ 0) > S^^{uj 0) 
for the LL parameters shown in Fig. |3.4| ). We obtain 

The theoretical shape of T'^ih^) plotted in FigJTsjat T = 250 mK > ^ J' 
is totally determined by the LL parameters (Fig. |3.4| ). Thus, similarly to the LL pa- 
rameters, it shows a strong asymmetry with respect to the center of the gapless phase 
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which is perfectly reproduced by the NMR measurements on BPCB (see Fig. |4.5|). 
The only free (scaling) parameter, A± = 0.041 T, is deduced from the fit of Eq. ( |4.3| ) 
to the experimental data f3T|. The fitted parameter agrees with that obtained from di- 
rect N^^ NMR determination | [3T| and rescale globally the theoretical curves without 
modifying its shape. For comparison, the T^^ obtained in the spin chain mapping 



approximation is also plotted in Fig. 4.5 As for other physical quantities, this de- 



scription fails to reproduce the non-symmetric shape. 



4.4 Properties of weakly coupled ladders 



As discussed in Sec. |2.1.2[ the interladder coupling J' induces a low temperature 



ordered phase (the 3D-ordered phase in Fig. 2.2 b). Using the mean field approxi- 



mation presented in Sec. 3.3 we characterize the ordering and compute the critical 



temperature and the order parameter related to this phase. 

4.4.1 3D order transition temperature 

In order to compute the critical temperature of the 3D transition, we follow Ref . 
and treat the staggered part of the mean field Hamiltonian Hmf ( |3.47| ) perturbatively 
using linear response. The instability of the resulting mean field transverse suscepti- 
bility for an order with momentum q = ir (staggered order), due to the 3D transition, 
appears at Tc when |99| 

X+-(<Z = 7r,c. = 0)|^^ = -^. (4.4) 

Where Xn" is the transverse correlation function of an isolated single ladder system 
defined in Eq. ( |3.33| ). As Tc ~ J' <C Jy, this correlation can be computed ana- 
lytically (see Eq. ( |3.36| )) using the LL low energy description of the isolated ladder 
(Eq. ( |3.24| )) in the gapless regime. Applying the condition ( |4.4| ) to the LL correla- 



tion ( [3.36| ) leads to the critical temperature 



^ 2iT \ 2u 




Introducing the computed LL parameters u, K and (see Fig. |3.4| ) in this expres- 
sion, we get the critical temperature | 31 1 as a function of the magnetic field. Only J' 
remains as a free parameter. Fitting our results for the transition temperature to the 
experimental data (Fig. |4.6| ) allows us to extract the mean field interladder coupling 
J(jp ^ 20 mK for the experimental compound BPCB. The asymmetry of the LL 
parameters induces a strong asymmetry of Tc with respect to the middle of the 3D 
phase which is in very good agreement with the experimentally observed asymmetry. 
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As the mean field approximation neglects the quantum fluctuations between the 
ladders, the critical temperature Tc is overestimated for a given J^p- In order to fully 
take the fluctuations into account a Quantum Monte Carlo (QMC) determination of 
Tc based on the same 3D lattice structure is performed in Ref. |38|. Let us note 
that QMC simulations of the coupled spin ladder Hamiltonian \2.l\ are possible 
since the 3D lattice structure, Fig. |2.6| is unfrustrated. Currently, this determination 
is only accessible for larger interladder couplings J' and shows |32| that the real 
critical temperature is well approximated by the mean field approximation, but with 
a rescaling of the real interladder coupling J' ^ 27 mK = J^p with a ^ 0.74. 
The rescaling factor a is essentially magnetic field independent [32J and similar to 
the values obtained for other quasi one-dimensional antiferromagnets |1QQ[|1QH . 




Figure 4.6: Magnetic field dependence of the transition temperature between the gapless 
regime and the 3D-ordered phase, Tc(/i^), is plotted in solid red line for the ladder LL 
parameters of BPCB shown in Fig. |3.4| (in dashed blue line for the LL parameters of the 
spin chain mapping). The NMR measurements from Ref. |31 ) are represented by black 
circles and the neutron diffraction measurements from Ref. 132) by green dots. [Taken 
from Ref. 13811 



4.4.2 Zero temperature 3D order parameter 



The staggered order parameter in the 3D-ordered phase, m^, can be analytically de- 
termined at zero temperature using the mean field appr oximation for the interladder 
coupling and the bosonization technique (see Sec. 3.3). As = y/A^{(zos{0{r))) 



ator 



9(r) 



in the bosonization description (Sec. 3.2) and the expectation value |96 | of the oper- 



is 



F{K) 



TTy/A^ricJ'ml 
2u 



(4.6) 
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for the sine-Gordon Hamiltonian Hsg ( |3.48| ) with 



8K 



F{K) 



S,K-1 



r(i-ak) 



\8K- 



1 J ^ \16K-2 J 



(4.7) 



we can extract 



(4.8) 



This can be evaluated in the 3D-ordered phase by introducing into (|4.8|) the LL pa- 



rameters li, and Ax from Fig. 3.4 Fig. 4.7 shows the order parameter versus the 



magnetic field determined analytically and numerically by DMRG (see Sec. 3.3.1 ). 
The two curves are almost indistinguishable and exhibit a strongly asymmetric camel- 
like shape | 3JJ with two maxima close to the critical fields. The asymmetry of the 
curve is again due to the presence of the additional triplet states. This asymmetry 
disappears in the spin chain mapping. 
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Figure 4.7: Magnetic field dependence of the transverse staggered magnetization per 
spin, rriaih^), at zero temperature in the 3D-ordered phase. Results of the analytical 
bosonization technique for the LL parameters of the compound BPCB shown in Fig. |3.4| 
and = 27 mK is represented by the dash-dotted red line (dashed blue line for the LL 
parameters of the spin chain mapping). The DMRG result for J' = 20 mK(= Jmf) is 
represented by black dots (as a comparison the bosonization result for = 20 mK is 
plotted in solid red line). Note, that these two curves are almost indistinguishable. A 
comparison of m^{h^) with the experimentally determined value using ND and NMR 



is shown in Fig. 4.12 [Taken from Ref. |^38J ] 
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4.5 Comparison with experimental results on BPCB 

Many experimental measurements related to the theoretical results presented at the 
beggining of this chapter have been performed on BPCB. In order to characterize this 
compound and understand its physical behavior, we discuss, in this section, several 
experiments performed on BPCB and compare these to the theoretical predictions. 

The longitudinal magnetization that can be measured very precisely by NMR at 
T = 40 mK (see Ref. |31 1) agrees remarkably well with the one computed using the 
weakly coupled ladder model (see Fig.|4.8|). In particular the linear growth close to 



the critical fields due to the small interladder coupling J' (Sec. 4.1 ) is highlighted (in 
the inset of Fig. |4.8| ). Nevertheless, the main shape of the magnetization is not very 
sensitive to the underlying model (see Fig. |4.1| and |4.8| ). Thus it cannot be used to 
distinguish between various models. However once the model is chosen e.g. a spin 
ladder, it can be used to fix precisely the parameters given the high accuracy of the 
experimental data. In particular, the position of the critical fields are very sensitive 



to the values of the intraladder couplings (Sec. 4.1 ). The couplings determined by 



the magnetization are ^ 12.6 K and Jm ^ 3.55 K. Note that the magnetization 



curve has been also measured by ND |32J (see Fig. |4.9| ) and agrees perfectly with 
the NMR experiment^ 

A more selective test to distinguish between various models is provided by the 
specific heat. This is due to the fact that the specific heat contains information on 
high energy excitations which are characteristic for the underlying model. As shown 



in Fig. |4.1Q| the experimental data from p3| are remarkably described, up to an 
accuracy of a few percent, by a simple Heisenberg ladder Hamiltonian with the pa- 
rameters extracted from the magnetization. In particular, not only the low tempera- 
ture behavior and the crossover from the LL regime to the quantum critical regime 



(Fig. |4.111 b) are covered by the ladder description, but also the higher maxima. This 
indicates that the ladder Hamiltonian is an adequate description of the compound and 
insure that no other large magnetic exchanges are forgotten. The small discrepancies 
between the specific heat data and the calculation which is essentially exact can have 
various sources. First of all, the substraction of the non-magnetic term in the ex- 
perimental data can account for some of the deviations. Furthermore the interladder 
coupling and the coupling anisotropics can induce slight changes in the behavior of 
the specific heat. 

The quality of the determination of the model and its intraladder parameters 
becomes more evident in the comparison of the NMR data for the relaxation rate 
T^^ with the theoretical results of the Luttinger liquid theory as shown in Fig. 4.5 



Only one adjustable parameter is left, namely the hyperfine coupling constant (see 



Sec. 4.3 ). This parameter allows one for a global rescaling of the theoretical curve, 



but not for a change of its shape which is totally determined by the LL parameters 



(Fig. 3.4 ). The agreement between the theory and the experimental data is very good 



: 2.17 for the experimental settings of ND measurements. 



47 



4. STATIC PROPERTIES AND NMR RELAXATION RATE 




Figure 4.8: Magnetic field dependence of the magnetization per Cu^^ ion in BPCB 
measured by NMR at T = 40 mK. The data are compared to the result of the 
DMRG calculation for a single ladder with the BPCB couplings (solid blue line) and 
for the spin chain mapping rescaled to fit with the single ladder critical fields (dashed 
green line), both at T = 0. Inset shows the critical linear dependence characteristic for 
weakly coupled ladders very close to /ic? (solid black line) and the determination of hc2 
using the assumption of a square root critical behavior for a single ladder (dashed pink 
line). See also Fig |4.1| for a comparison of the theoretical computations. [Taken from 
Ref. ^31J] 
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Figure 4.9: Magnetic field dependence of the square of the magnetization per rung in 
BPCB, m^^(/i^), measured by ND at T = 54 mK (blue circles). The solid red line 
represents the theoretical prediction computed by DMRG with the BPCB couplings. 
[Taken from Ref. |j32J] 
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Figure 4.10: Temperature dependence of the specific heat measurements c(T) on the 
compound BPCB from Ref. 1 33 1 (dots) and the T-DMRG calculations (solid lines) in the 
spin liquid regime at (a.l) /i^ = T, (a.2) /i^ = 3 T, (a.3) /i^ = 5 T and in the gapless 
LL regime at (b.l) /i^ = 9 T, (b.2) = 10 T, (b.3) /i^ = 11 T. See also Fig.|^for a 
comparison of the theoretical computations and Fig. 2.3 for a false color picture of c/T 



in a full range of temperature and magnetic field. [Data taken from Ref. | 38 1] 
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Figure 4.11: Magnetocaloric effect measured on BPCB. (a) Heat-flow to and from 
the sample as a function of magnetic field divided by temperature, (6Q/Sh^)/T = 
— {dM^/dT)\hz. (b) Crossover temperature Tll of the LL to the quantum critical 
regime versus the applied magnetic field . White circles (black circles) denote the 
phase boundary derived from the {dm^ / dT)\h^ = criterium computed from BPCB 
measurements shown in (a) ({dc/dT) | = criterium computed from BPCB measure- 
ments shown in Figs. |4.10"| and |2.3| l. As a comparison the blue circles and red crosses 
show the T-DMRG computations of this crossover presented in Fig. |4.3| for both criteria 
respectively. [Taken from Ref. |^33J] 
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over the whole range of the magnetic field and only small deviations can be seen. 
As discussed further, two other correlations included in the LL description of the 
low temperature 3D order and its critical temperature are tested with the same LL 
parameters. This compound thus allows us to quantitatively test the Luttinger liquid 
universality class. Even though the Luttinger liquid description is restricted to low 
energies, in BPCB its range of validity is rather large. Indeed at high energy, its 
breakdown is approximately signaled by the first peak of the specific heat 1 33 1 (see 



Sec. 4.2.2) or the cancellation of the magnetocaloric effect (equivalent to the first 



extrema of the magnetization versus the temperature discussed in Sec. |4.2.1| ). Here 
the experimentally determined crossover is located about T ^ 1.5 K at midpoint 
between hd and hc2, and agrees totally with its numerically computed value (see 



Fig. 4.11 b, respectively). Given the low ordering temperature which has a max- 
imum at about T ~ 100 mK this leaves a rather large Luttinger regime for this 
compound. 

Finally, deviations from the simple ladder Hamiltonian can be present. Small 
anisotropy of the couplings can exist and indeed are necessary to interpret recent ESR 
experiments | [5Q| . Other terms such as longer range exchanges or Dzyaloshinskii- 
Moryia (DM) terms might occur along the legs even if the latter is forbidden by 
symmetry along the dominant rung coupling. Clearly all these deviations from the 
Heisenberg model cannot be larger than a few percents. They will not lead to any 



sizeable deviation for the Luttinger parameters (Fig. |3.4| ) in the one dimensional 
regime. Close to the critical points they can, however, play a more important role. 
It would thus be interesting in subsequent studies to refine the model to take such 
deviations into account. 

Taking now the coupling between ladders into account, one can induce a tran- 
sition to a three-dimensional ordered phase. The transition temperature is shown in 
Fig. |4.6[ Experimentally it is determined by NMR | [3T| and neutron diffraction mea- 



surements |32|. Theoretically the ladders are described by Luttinger liquid theory 
and their interladder coupling is treated in a mean field approximation (Sees. |3.2| 



and |3.3| ). As shown in Fig. |4.6| the Luttinger liquid theory provides a remarkable 
description of the transition to the transverse antiferromagnetic order at low temper- 
atures. The shape of Tdh^) is almost perfectly reproduced, in agreement with both 
the NMR | [3T| and the ND data p2| . The comparison with the experiments deter- 
mines the interladder coupling J' , the only adjustable parameter. The simple mean 



field approximation would give a value of J' ~ 20 mK. As discussed in Sec. |4.4.1 
mean field tends to underestimate the coupling and it should be corrected by an es- 
sentially field independent factor. Taking this into account we obtain a coupling of 
the order of J' = 27 mK. 

The order parameter in the antiferromagnetic phase can also be observed by ex- 



periments. As discussed in Sec. 4.4.1 it shows a very interesting shape. At a pure 



experimental level neutron diffraction and NMR have some discrepancies as shown 



in Fig. 4.12 These discrepancies can be attributed to the different temperatures at 



50 



4.5 Comparison with experimental results on BPCB 
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Figure 4.12: Magnetic field dependence of the transverse staggered magnetization per 
spin, rriaih^), in the 3D-ordered phase. Comparison between NMR measurements 
(black circles) done at T = 40 mK from Ref. | 31 1 and scaled to the theoretical results 
for J' = 27 mK (dash-dotted red line), neutron diffraction measurements on an absolute 
scale from Ref. |32| at T = 54 mK (T = 75 mK) (red crosses (black dots)). Recent 
neutron diffraction measurements as a function of temperature suggest that the data of 
Ref. 1 32 1 was taken at temperatures approximately 10 mK higher than the nominal in- 
dicated temperature. See also Fig |4.7| for a comparison of the theoretical computations. 
[Taken from Ref . | 38 |] 

which the data have been taken, and a probable underestimation of the temperature 
in the neutron diffraction experiments |32|. Indeed the order parameter close to the 
critical magnetic field hc2 is very sensitive to temperature, since the transition tem- 
perature drops steeply in this regime. Note that although the NMR allows clearly 
for a more precise measurement of the transverse staggered magnetization it cannot 
give its absolute value. Thus the amplitude of the order parameter is fixed from the 
neutron diffraction measurement. Even though a good agreement between the theo- 
retical results and the experimental results is obtained, several questions concerning 
the deviations remain to be addressed. 

First, the theoretical curve does not fully follow the shape of the experimental 
data. Particularly at high fields the experimental data shows a stronger decrease. A 
simple explanation for this effect most likely comes from the fact that the calculation 
is performed at zero temperature, while the measurement is done at 40 mK. This is 
not a negligible temperature with respect to Tc, in particular at magnetic fields close 
to hc2- Extrapolation of the experimental data to zero temperature | [3T| improves 
the agreement. Nevertheless, for a detailed comparison either lower temperature 
measurements or a calculation of the transverse staggered magnetization at finite 
temperature would be required. From a theoretical point of view, including the 3D 
coupling in the model ( |2.1| ) make the finite temperature computations difficult to 
perform (analytically and numerically). Such an investigation thus clearly require 
the development of more powerful techniques. 
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The second question comes from the amphtude of the staggered magnetization. 
Indeed the experimental data seem to be slightly above the theoretical curve, even 
if one uses the value J' = 21 mK for the interladder coupling. This is surprising 
since one would expect that going beyond the mean field approximation could only 
reduce the order parameter. Naively, one would thus need a larger coupling, perhaps 
of the order of J' ^ 60 — 80 mK to explain the amplitude of the order parameter. 
This is a much larger value than the one extracted from the comparison of Tc. How 
to reconcile these two values remains open. The various anisotropics and additional 
small perturbations in the ladder Hamiltonian could resolve part of this discrepancy. 
However, it seems unlikely that they result in a correction of J' by a factor of about 
2-3. Another origin might be the presence of some level of frustration present in the 
interladder coupling. Clearly more experimental and theoretical studies are needed 
on that point. 
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CHAPTER 5 



Dynamical correlations of a spin ladder 



In this chapter, we focus on the zero temperature spectral functions of a single spin- 



1/2 ladder computed with t-DMRG (see Sec. 3.1.3 ). As we will see in the following 
these dynamical quantities are direct probes of the excitations in the system and are 
experimentally accessible through INS measurements. 

After an introduction of the relevant spectral functions related to spin- 1/2 lad- 
ders, we discuss the possible rung excitations created by the spin operators. Next, 
we analyze in detail the computed spectra for the parameters of the compound BPCB 



(see Sec. |2.2| ) separately in the gapped spin liquid and the gapless regime. In addi- 
tion, these spectra are compared to analytical results when such results exist. In 
particular, we check the agreement with the LL description at low energy and use 
a strong coupling expansion (appendix [A|) to qualitatively characterize the different 
excitations occurring. Then we generalize our analysis for different coupling ratios 
7 from the weak (7 00) to strong coupling (7 ^ 0), and briefly discuss the in- 
fluence of the weak interladder coupling on the excitations of the system. Finally, in 



Sec. 5.6 we compare the low energy part of the computed spectra to the INS mea- 
surements on BPCB and provide a quantitative prediction for the high energy part of 
the INS spectra. Furthermore we give a short overview of the ND experiments for 
the measurement of static quantities. 
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5.1 Zero temperature spectral functions 

As discussed in Sec. |3.2.3[ in a ladder system different types of correlations are 
possible. We focus here on the quantities 

/oo 
(5.1) 

where 6'^^^^ = Sf"^-^ ± ^[^2 the symmetric (+) and antisymmetric (— ) operators 
with rung momentunj^ g'^/ = parity in the rung direction P = +1,-1, 

respectively. The type of the correlation is denoted by a, /3 = 2:, +, — . The time 
evolution {t) = e^^^S'^^^^e"*^^ is with respect to the Hamiltonian H ( |2.3| ) of a 
single spin ladder. 

These correlations are lattice versions of the spectral functions introduced in 



Sec. 3.2.3 with qy and q the momenta in the rung and along the leg lying in the 
first Brillouin zone. The latter is given in reciprocal lattice units a~^. 

Using the reflection and translation invariance of an infinite size system (L 
00), we can rewrite the considered correlations ( |5.1| ) with S^^ = 5*^ in a spectral 
decomposition (at zero temperature), i.e. 

= ^ ^ \{X\Sl{qm\'S{uj + Eo- Ex) (5.2) 

A 

where |0) denotes the ground state of H with energy Eq, S^^{q) = 

the sum over all eigenstates |A) of i^, and Ex their eigenenergy. The form of 
Eq. \5.2) clearly shows that Sq^{q^uj) is non-zero if the operator S^^ can create 
an excitation |A) of energy Eq -\- uj and momentum q from the ground state. The 
correlations S^^ are then direct probes of the excitations | A) in the system. 

Since the experimentally relevant case (compound BPCB) corresponds to a rel- 
atively strong cou pling s ituation (7 <C 1, Eq. |2.6| ), we use the decoupled bond limit 
to represent the expected excitations on a single rung 



introduced in Sec. 



2.1.1 



\t^), \t ) or l^). In table 5.1 we summarize the rung excitations created by the 
operators S^^ and their properties. For example the operator applied on the sin- 
glet \s) excites the triplet \t^). Typically the rung parity P is changed by applying 
an operator with rung momentum qy = tt and the z-magnetization is modified by 
AM^ = ±1 by applying the operators , respectively. 



5.2 Excitations in the spin liquid 

Using the decoupled bond limit in the spin liquid phase, the excitations in the system 
can be pictured as the excitation of rung singlets to rung triplets. At zero magnetic 

^The rung momentum qy is a good quantum number. 
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Table 5.1: Rung excitations created by the symmetric and antisymmetric operators in 
the decoupled bond Umit. The elements of the first column represents the initial rung 
states on which the rung operators written in the first line apply. The effect of of these 
operators on the parity P and the magnetization of the system is also summarized 
in the two last lines. 



field = 0, the system is spin rotational symmetric and the different triplet ex- 
citations have the same energy ~ Jj_. It has been seen previously that in the spin 
liquid both single triplet excitations and two-triplet excitations play an important role 
| T4f|T8| . We discuss these excitations in the following focusing on the ones that can 
be created by the symmetric Sq^ = 2S^^ and the antisymmetric S^^ = 2S^^ cor- 
relations (see Fig. |5.1| ) for the BPCB parameters ( |2.15| ), where a = x^y^ z . Note that 
these correlations are independent of the direction due to the spin rotation symmetry. 



5.2.1 Single triplet excitation 

At zero magnetic field = 0, the system is in a global spin singlet state (S* = 0) p|. 
The Qy = TT correlation couples this ground state to states with an odd number of 
triplet excitations with rung parity P = — 1 and total spin S = 1^ = ±1,0 (see 
table |5.1| ). Nevertheless, only single triplet excitations are numerically resolved as 



shown in Fig. 5.1b. Their spectral weight is concentrated in a very sharp peak whose 
dispersion relation, uJt{q), can be approximated using a strong coupling expansion in 
7 similar to the one described in appendix |a| Up to thirc^ order y_5j| the dispersion 
is given by 

l+7cos(7+^(3-cos2(7) + ^(3-2cos(7-2cos2(7 + cos3(7) + 0(7'^). 
J±_ 4 8 

(5.3) 

At the first order, it is simply a cosine dispersion [14J , i.e. uJt{q)/J± ~ 1 + 7 cos q. 
The lowest energy single triplet excitation has a momentum q = ir and corresponds 



■ Note that an expression up to seventh order in 7 has been determined in Ref. |102| 
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71/2 71 371/2 271 7r/2 % 37i/2 27i 

q(a-^) q(a-^) 

Figure 5.1: Momentum-energy dependent correlation functions Sq^^q^uo) with a — 
x,y, z at = 0. Numerical results are shown with the color coding, (a) Symmetric 
part So^{q,(jj). The dashed (black) line marks the (q^oj) position of the two-triplet 
bound state with dispersion relation uJtB{q), Eq. ( |5.4| ), in its existence interval ( |5.5| ). 
The dash-dotted (white) lines correspond to the boundaries of the two-triplet continuum, 
(b) Antisymmetric part S!^"^ {q, cj). The dashed (black) line corresponds to the predicted 
dispersion relation (jOt{q) of a single triplet excitation, Eq. \53\ . [Taken from Ref. |38|] 



2.1.1 



and 4.1 h. In 



to the spin liquid ga^ or the first critical field hd (see Sees. 
Fig. 5.1 b we compare the numerical results for the BPCB parameters ( |2.15| ) to the 
expression ( |5.3| ). The strong coupling expansion describes very well the position of 
the excitations found numerically. The comparison with the known solutions serves 
as a check of the quality of our numerical results. 



5.2.2 Two-triplet excitations 



The structure of the Qy = correlation is more complex (Fig. 5.1 a). Due to the rung 
parity P = 1 of the operators Sq, the excitations correspond to an even number of 
triplet excitations with total spin 6* = 1, = ±1,0 excited from rung triplets 
already present in the ground state. We focus here on the two-triplet excitations that 



can be resolved numerically (Fig. 5.1 a). These can be divided into a broad contin- 
uum and a very sharp triplet (S = 1) bound state of a pair of rung triplets. Since 
these excitations stem from the coupling to triplets already present in the ground 



state (Fig. |4.2| ), their amplitude for the considered BPCB parameters ( [2.15| ) is con- 



siderably smaller than the weight of the single triplet excitations p8| . 

The dispersion relation of the bound states, CjOtsiQ), has been calculated using a 



^Note that the spin gap has been determined up to 13th order in 7 (121. 
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linked cluster series expansion [16} up to third order in 7 

^tB{q) _ 2 + 2(-3 - 2cos(7) + ^(11 - 2 cos g - 4 cos 2^) 



J_L 2 



(17 + 9 COS (7 - 8 cos 2q - 5 cos 3q) + 0{j^). (5.4) 

16 



The first terms of the expansion have an inverse cosine form and the bound state only 
exists in an interval 

,.<,<2.-,. with ^c=y-^-^ + 0(7^) (5.5) 



around q = n (cf. Ref. |T6]T7|). The numerical results for the BPCB parame- 



ters ( |2.15| ) agree very well with the analytic form of the dispersion (Fig. |5.1[ a). 
The upper and lower limits of the continuum can be determined by considering the 
boundary of the two non-interacting triplet continuum consisting of two excitations 
with the single triplet dispersion ( |5.3| ). They agree very well with the results found 
numerically (Fig . |5.1[ a) . 



5.3 Excitations in the gapless regime 

A small applied magnetic field (h^ < hd), at first order, smoothly translates the 



excitations shown in Fig. 5.1 by an energy —h^M^, with the magnetization of 
each excitation, due to the Zeeman effect. However, if the magnetic field exceeds 
hci, the system enters into the gapless regime and the structure of the excitations 
spectrum changes drastically. A continuum of excitations at low energy arises. For 
small values of 7 most features of this low energy continuum are qualitatively well 
described by considering the lowest two modes of the ladder only. Beside the low 
energy continuum, a complex structure of high energy excitations exist. Contrarily 
to the low energy sector, this structure crucially depends on the high energy triplet 
modes. In the following, we give a simple picture for the excitations starting from 
the decoupled bond limit. 

5.3.1 Characterization of the excitations in the decoupled bond 
limit 

The evolution of the spectra for the BPCB parameters with increasing magnetic field 
are presented in Fig. 5.2 for S^^, in Fig. 5.3 for S^~, and in Fig. 5.4 for S~^. Three 



different classes of excitations occur: 

(i) a continuum of excitations at low energy for Sq^ and S^^ 
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Figure 5.2: Momentum-energy dependent zz-correlation function (1) at — 0.25 
{h^ = 3.153 J||), (2) at = 0.5 (/i^ = 4.194 Jy), and (3) at = 0.75 
{h^ — 5.192 J||). (a) Symmetric part S^^i^q^uS) without Bragg peak at ^ = 0. The 
dashed black Unes correspond to the location of the slow divergences at the lower edge 
of the continuum predicted by the LL theory. The dashed white curve corresponds to 
the predicted two-triplet bound state location, (b) Antisymmetric part Sl^{q^ uj). The 
dashed black lines correspond to the position of the high energy divergences or cusps 
predicted by the approximate mapping onto the t-J model. The vertical white dash- 
dotted lines mark the momenta of the minimum energy of the high energy continuum 
and the black cross is the energy of its lower edge f22 1 at ^ = tt. [Taken from Ref. f38l] 
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Figure 5.3: Momentum-energy dependent H correlation function (1) at = 0.25 

(h^ = 3.153 J||), (2) at = 0.5 (/i^ = 4.194 J\\), and (3) at = 0.75 (/i^ = 
5.192 J||). (a) Symmetric part Sq {q, cj). The vertical dash-dotted white lines mark the 
momenta of the minimum energy of the high energy continuum and the horizontal ones 
the frequency of its lower edge |22 | at ^ = 0, 27r. The dashed white lines correspond 
to the position of the high energy divergences or cusps predicted by the approximate 
mapping onto the t-J model. The dotted white curve corresponds to the predicted two- 
triplet bound state location at cj 3hz which is hardly visible, (b) Antisymmetric part 
Si~ {q, (jj). The dashed and dash-dotted (dotted) white lines correspond to the location 
of the strong divergences (cusps) at the lower edge of the continuum predicted by the LL 
theory. [Taken from Ref. 138)] 
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Figure 5.4: Momentum-energy dependent — h-correlation function (1) at = 0.25 
(h^ = 3.153 J||), (2) at = 0.5 (h^ = 4.194 Jy), and (3) at = 0.75 
(h^ = 5.192 J||). (a) Symmetric part SQ~^{q,uj). The vertical dash-dotted white Unes 
correspond to the momenta at which the minimum energy of the high energy contin- 
uum occurs and the horizontal line to the frequency of its lower edge | 22 1 at ^ = 0, 27r. 
The dashed black curve corresponds to the predicted two-triplet bound state location, 
(b) Antisymmetric part S~~^{q,(jj). The dashed and dash-dotted (dotted) white lines 
correspond to the location of the strong divergences (cusps) at the lower edge of the 
continuum predicted by the LL theory. [Taken from Ref. ||38j|] 
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(ii) single triplet excitations at higher energy with a clear substructure for S^^ , 

5'+-,and5+- 

(iii) excitations at higher energy for Sq^ and Sq~ and 5*0"^ stemming from two- 
triplet excitations which have their main weight around it. 

In the following we summarize some of the characteristic features of these excita- 



tions, before we study them in more detail in Sees. 5.3.2 and 5.3.3 



(i) The continuum at low energy which does not exist in the spin liquid is a char- 
acteristic signature of the gapless regime. It stems from excitations within the 
low energy band which corresponds to the \s) and states in the decoupled 
bond limit (cf. Fig.jZ^ja and table [51]): 

Sq^ : excitations within the triplet mode 
S^^ : excitations between the singlet \s) and the triplet mode. 

This continuum is smoothly connected to the spin liquid spectrum in the case 



of . It originates from the single triplet branch (Sec. 5.2.1) when 
the latter reaches the ground state energy due to the Zeeman effect. Since two 
modes play the main role in the description of these low energy features, many 
of them can already be explained qualitatively by the spin chain mapping. The 
excitations in the chain have been studied previously using a Bethe ansatz de- 
scription and exact diagonalization calculations in Ref. 1 103|. More recently 
they were computed in Ref. 1 104] to a high accuracy due to recent progress in 
the Bethe ansatz method. In particular, the boundary of the spectrum at very 
low energy is well described by this approach, since the LL velocity determin- 
ing it is hardly influenced by the higher modes (cf. Fig. |3.4| ). However, a more 
quantitative description requires to take into account the higher modes of the 



system as well. In Sec. 5.3.2 we compare in detail our results with the LL 
theory and the spin chain mapping pointing out their corresponding ranges of 
validity. 

(ii) The single high energy triplet excitations form a continuum with a clear sub- 
structure. In the decoupled bond limit, these excitations correspond to 



5- 

s+- 



Single triplet excitations \t^) at energy 
Single triplet excitations \t^) at energy 
Single triplet excitations \t~) at energy 2h^ . 



Many of the features of these continua can be understood by mapping the 
problem onto a mobile hole in a spin chain, as pointed out first in Ref. | |1Q5| . 
We detail in Sec. 5.3.3.2 and appendix [A| this mapping. It opens the possibility 



to investigate the behavior of a single hole in a t- J like model using experiments 
in pure spin ladder compounds. 
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(iii) The high energy continuum, which has almost no weight close to the Brillouin 
zone boundary = 0, 27r), is related to two-triplet excitations of the spin 
liquid (Sec. |5.2.2| ). They are generated from high energy triplet components of 
the ground state. Their weight therefore vanishes for 7 ^ and the excitations 
correspond to 

\t^)\t^)) at energy - 
\t-)\t^)) at energy - 2h^ 
|t-)|t°))atenergy-3/i^ 

5.3.2 Low energy continuum 

In this section we concentrate on the low energy excitations of type (i) discussing 
first their support and then comparing their spectral weight to the LL prediction. 



Sq^ : Two-triplet excitations - 
: Two-triplet excitations -^{\t^)\t~) - 
: Two-triplet excitations -j^{\t^)\t~) - 



5.3.2.1 Support of the low energy excitations 



The position of the soft modes in the low energy continuum can be directly obtained 
from the bosonization representation |[3][2TJ[22] (see Sec. |3.2.3) . They can also be 



understood in a simple picture which we outline in the following. The distribution of 
the rung state population in the ground state depends on the magnetic field (see 



Fig. |4.2| ). Taking a fermionic point of view, the magnetic field acts as a chemical po- 
tential that fixes the occupation of the singlet and triplet rung states. Increasing the 
magnetic field reduces the number of singlets, whereas at the same time the number 



of triplets increases (see sketch in Fig. 5.5 ). The Fermi level lies at the momenta q 



7rm^ , TT ( 2 — ) for the singlet states and at the momenta g = 7r(l — m^),7r(l + m^) 
for the triplet states. In this picture the soft modes correspond to excitations at the 
Fermi levels. For transitions ^ the transferred momenta of these zero 
energy excitations are g = 0, 27rm^, 27r(l — m^). In contrast, the interspecies tran- 
sitions ^ \s) allow the transfer of g = 7r(l — 2m^), tt, 7r(l + 2m^). Therefore, 
the positions of the soft modes in the longitudinal correlation Sq^ which correspond 
to transitions within the triplet states shift from the boundaries of the Brillouin zone 



inwards towards q = n when increases (Fig. 5.2 a). In contrast, the positions of 



the soft modes in the transverse correlations S^^ which induce transitions between 



the singlets and the triplets move with increasing magnetic field outwards (Figs. 5.3 b 
and[54lb). 

The top of these low energy continua are reached when the excitations reach the 
boundaries of the energy band. In particular, the maximum of the higher boundary 
lies at the momentum q = ir which is easily understood within the simple picture 



drawn above (cf. Fig. 5.5). A more detailed description of different parts of these 



low energy continua is given in Ref. | |103| |. 
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(a) (b) 




Figure 5.5: Fermionic picture for the effect of the magnetic field from : FilUng of (a) the 
singlet band |s), (b) the triplet band |t^) in the gapless phase for a given magnetization 
m^ [Taken from Ref . |38)] 



Let us compare the above findings with the predictions of the LL theory for the 
dynamical correlations |[3j|2T]|22|. Details on the LL description of the correla- 
tions are given in Sec. |3.2.3l The LL theory predicts a linear momentum-frequency 
dependence of the lower continuum edges with a slope given by the LL velocity 
±.u (Fig. |3.4|). The position of the soft modes are given by the ones outlined in 



Sec. 3.2.3.2| (see Fig. 3.5). The predicted support at low energy agrees very well 



with the numerical results (Fig. |5.2[ a, |5.3[ b, and |5.4[ b). Of course when one reaches 
energies of order Jy in the spectra one cannot rely on the LL theory anymore. This 
is true in particular for the upper limit of the spectra. 



5.3.2.2 Spectral weight of the low energy excitations 

Let us now focus on the distribution of the spectral weight in the low energy con- 
tinuum. In particular, we compare our numerical findings to the Luttinger liquid 
description. Qualitatively, the LL theory predictions for the low energy spectra are 
well reproduced by the DMRG computations. 

The Luttinger liquid predicts typically an algebraic behavior of the correlations 
at the low energy boundaries which can be a divergence or a cusp. 

Sq^ : The Luttinger liquid predicts peaks at the g = 0, 27r branches and a slow diver- 
gence at the lower edge of the incommensurate branches q = 27rm^, 27r(l — 



m^) (with exponent 1 — K ^ 0.2 <C 1 (Fig. 3.6 b)). In the numerical re 



suits (Fig. 5.2 a) a slight increase of the weight towards the lower edge of the 



incommensurate branches can be seen. 
: A strong divergence at the lower edge of the q = tt branch (with exponent 



1 — 1/4K ^ 3/4 (Fig. 3.6 a)) is obtained within the Luttinger liquid de- 



scription. This is in good qualitative agreement with the strong increase of the 



spectral weight observed in the numerical data (Fig. 5.3 ). A more interesting 
behavior is found close to the momenta g = 7r(l ± 2m^) in the incommen- 
surate branches. Here a strong divergence is predicted for momenta higher 
(lower) than the soft mode q = 7r{l — 2m^) (q = 7r(l + 2m^)) with expo- 
nent 1 — T]- ^3/4::^0 (Fig.|3.6[a). In contrast for momenta lower (higher) 
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than the soft mode g = 7r(l — 2m^) (q = 7r(l + 2m^)) a cusp with expo- 
nent 1 — 7^+^— 5/4<C0 (Fig. 3.6 c) is expected. In the numerical results 



(Fig. |5.3[ b) this very different behavior below and above the soft modes is 
evident. The divergence and cusp correspond to a large and invisible weight, 
respectively. 



-+ . 



The same behavior as for 5*+ replacing 
Fig.[54lb. 



-m^ can be observed in 




Figure 5.6: Cuts at fixed momentum q = tt and magnetization = 0.5 of the low 
energy spectrum (a) So^(q = tt^lj), and (b) S~~^{q = tt^u). The (red) circles and 
the (black) squares are the numerical results for the ladder and its spin chain mapping, 
respectively The dashed lines correspond to the LL predictions and the solid lines are 
the latter convolved with the same Gaussian filter than the numerical data. The DMRG 
frequency numerical limitation is of the order of the peak broadening of width 6uj ~ 
0.1 J|| (see Sec. [3XTT| . [Taken from Ref [38 1] 



To compare quantitatively the predictions of the LL to the numerical results we 
show in Fig. |5.6| different cuts of the correlations at fixed momentum q = n and 
magnetization = 0.5 for the ladder and the spin chain mapping. These plots 
show the DMRG results, the LL description, and the latter convolved with the Gaus- 
sian filter. The filter had been used in the numerical data to avoid effects due to 



the finite time-interval simulated (see Sec. 3.1.3.1 ). Note that the amplitude of the 
LL results are inferred from the static correlation functions, such that the LL curve 
is fully determined and no fitting parameter is left. Therefore, the convolved LL 
results can directly be compared to the numerical results. Even though the actual 
numerical resolution might not be good enough to resolve the behavior close to the 
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divergences (cusps), interesting information as the arising differences between the 
spin chain mapping and the full ladder calculations can already be extracted. In 
Fig. 5.6 a, we show a cut through the correlation SQ^{q = 7r,cj). The convolved 
LL and the numerical results compare very well. The difference between the real 
ladder calculations and the spin chain mapping that neglect the effects of the higher 
triplet states \t~), \t^) is obvious. From the LL description point of view, the shift of 
the spin chain correlation compared to the real ladder curve comes mainly from the 
prefactor Az and the algebraic exponent which are clearly modified by the effects of 
the high energy triplets (see Fig. |3.4| ). 

For the transverse correlations, the LL theory predicts a strong divergence (with 
an exponent 1 — 1 / ^ 3/4 at the lower boundary of the continuum branch 
at g = TT. A cut through the low energy continuum S~~^{q = ir^u) is shown in 
Fig. 5.6 b. The convolved Luttinger liquid reproduces well the numerical results. 

From the comparisons in Fig. |5.6| we show that the numerical computations al- 
lowing one to access the middle and high energy excitations and the analytical LL 
theory describing the low energy physics are complementary. Due to their large over- 
lap, the combination of the two approaches provide a quantitative description of the 
correlations in a full range of energy. Nevertheless, a refined numerical investiga- 
tion of the effects close to the continuum boundaries would be required in order to 
investigate the non LL edge singularities recently pointed out in Ref. | |lQ6f|lQ8| . 



5.3.3 High energy excitations 

Before looking in detail at the two kinds of high energy excitations presented in 
Sec. |5.3.l] we compare briefly our computed high energy spectra with the weak cou- 
pling description (7 1). 



5.3.3.1 Weak coupling description of the high energy excitations 

In the weak coupling limit, information on the spectrum can be extracted from the 
bosonization description | 3||2Tp2l. In particular one expects a power law singularity 
at the lower edge continuum with a minimal position at g = 7r(l ± m^) (for S^^), 
q = 7rm^, 7r(2 — m^) (for Sq^) and an energy at momentum q = n (for S^^), 
q = (for S^^). Except for 5*0"^ in which the spectral weight is too low for a good 
visualization, our computed spectra reproduces well the predictions for the minimal 
positions even though the coupling strength considered is not in the weak coupling 
limit (cf. Figs.|5.2[b,|5.3[a and|5.4|a). 



5.3.3.2 High energy single triplet excitations 

The high energy single triplet continua originate from the transition of the low energy 
rung states \s) and to the high energy triplets \t^) and \t~). The excitations 
coming from the singlets \s) (in S^^ and S:^~) are already present in the spin liquid 
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phase (cf. Sec. |5.2.1| ) in which they have the shape of a sharp peak centered on the 



triplet dispersion. The transition between the gapped spin liquid and the gapless 
regime is smooth and consists in a splitting and a broadening of the triplet branch 
that generates a broad continuum of new excitations. Contrarily to the latter the 
excitations coming from the low energy triplets (in Sq~) are not present in the 
spin liquid phase. The corresponding spectral weight appears when > hd . 

An interpretation of the complex structure of these high energy continua can be 
obtained in terms of itinerant quantum chains. Using a strong coupling expansion 
of the Hamiltonian ( |2.3| ) (appendix |A|) one can map the high energy single triplet 
excitations \t^) to a single hole in a system populated by t wo type s of particles with 



2.1.1 



pseudo spin |f ) = \l) = \s) (with the notation of Sec. ! 

In this picture the effective Hamiltonian of the J± energy sector is approxi- 
mately equivalent to the half filled anisotropic ID t-J model with one hole (see ap- 
pendix |A3]2]). The effective Hamiltonian is given by 

Ht.j = Hxxz + i^t + i^s-h + e. (5.6) 

where e = ( + /i^)/2 is an energy shift, Hxxz is the XXZ spin- 1/2 chain Hamil- 
tonian ( [2.11| ) and = J||/2^^ ai^l a^i-\-i a + ^-^O usual hopping term. 
Here cj ^ (ci^a) is the creation (annihilation) operator of a fermion with pseudo spin 
(T = t, I at the site Note that although we are dealing here with spin states, it is 
possible to faithfully represent the three states of each site's Hilbert space (\s), 
\t^)) using a fermionic representation. 

In addition to the usual terms of the t-J model, a nearest neighbor interaction 
term between one of the spins and the hole arises 



J\\ ^--v r 



(5.7) 



Here ni^h^^ the density operator of the hole on the site 1. In this language the spectral 
weight of S^^ and Sq~ corresponding to the single high energy triplet excitations 
is equivalent to the single particle spectral functions of the spin up and spin down 
particle, respectively: 

5'f (X (etc-) with hole of type \s) ^ \t^) 

5'+" (X (etc-) with hole of type \t^)^\t^). ^^'^^ 

Here {cic,){q^uj) = Ea |(A|c,,,|0)p(^(a; + - ^a). 

For the standard t-J model (for SU(2) invariant XXX spin chain background and 
without the anisotropic term i^s-h in Eq. ( |5.6| )), these spectral functions have been 
studied in Refs. |1Q9||11Q| . The presence of singularities of the form 

{clc,){q,uj) oc [c-c.to(g-g.)]2^^(«)-i (5.9) 
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were found. Here cj^o (q) is the \t^) triplet dispersion relation, Qjy the spinon momen- 
tum at the Fermi level and the algebraic decay exponent at the singularity. This 
exponent is not known in our case and depends on the magnetization and the 
momentum q. Eq. ( |5.9| ) generates a peak or a cusp at the energy uj = uoto {q — Qu)- 
The spinon momentum Qj^ depends on the type of the rung state before excitation 
(u = 5, For an excitation created from a singlet state one obt ains Qs = ±7rm^ 
(for S^^) and from the triplet state qt+ = 7r(l ±m^) (for S^~) (Fig.jsis]). At = 0, 
a series expansion of cj^o (q) can be performed (cf. ujt{q) in Eq. ( |5.3| )). To extend this 
expression into the gapless phase (hd < < hc2), we approximate uoto{q) by 
shifting the value ujt{q) at /i^ = by the Zeeman shift, i.e. 



(5.10) 



Here we used the shift of the ground state energy per rung AEo{h^) = EQ{h^) — 
Eo{0). AEo was determined by DMRG calculations (Fig.|5^for the BPCB param- 
eters). The resulting momentum- frequency positions uj = uJto{q — qj^) of the high 



energy singularities (cusps or divergencies) are plotted on the spectrum Fig. 5.2 b 
and Fig. |5.3| a. They agree remarkably well with the shape of the computed spectra, 
in particular, for small magnetic field |^ Neglecting the additional interaction term 
i^s-h, the t-J model Hamiltonian would lead to a symmetry of these excitations with 
respect to half magnetization. However, in the numerical spectra the effect of the in- 
teraction shows up in a clear asymmetry of these excitations (compare Figs. |5.2| l.b 
and 5.3 3. a). In particular, in the ^q"" correlation some of the weight is seemingly 
detaching and pushed towards the upper boundary of the continuum (Fig. |5.3| 3.a) 
for large magnetization. 
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Figure 5.7: Shift of the ground state energy per rung versus the magnetization 
A^o(m^) = ^o(O) - Eo{m^). [Taken from Ref . |38)] 

A similar mapping can be performed for the single \t~) excitation. In contrast to 
the Jj_ sector, in the 2 sector not only the \t~) excitation occurs, but the effective 
Hamiltonian mixes also \t^) triplets into the description. Therefore, the description 

^For the correlations S^^ and Sq , some of these singularities correspond to the lower edge descrip- 
tion in Ref. \22j and discussed in Sec. |5.3.3TT] 
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by a single hole in a spin- 1/2 chain breaks down and more local degrees of freedom 
are required. This results in a more complex structure as seen in Fig. 5.3 l.b. Pre- 
viously high-energy excitations in dimerized antiferromagnets have been described 
by a rather general mapping to an X-ray edge singularity problem |lllf[113| . It is 



interesting though that in the present setup these excitations can be understood as t-J 
hole spectral functions, which display a much richer structure than anticipated. 



5.3.3.3 High energy two-triplet excitations 

The two-triplet continua and bound states already discussed in the spin liquid phase 
(cf. Sec. |5.2.2| ) are still visible in the gapless regime in the symmetric correlations 
(Sq^ and S^^). At low magnetic field the location of their maximal spectral weight 
can be approximated by the expression of the bound state dispersion at zero field, 
(jOtsiQ) in Eq. ( |5.4| ), shifted by the Zeeman energ}j^ The two-triplet excitation lo- 
cation obtained in this way agrees to a good extent with the location found in the 
numerical calculations (cf. Figs. |5.2[ l.a, |5.3[ l.a and |5.4[ l.a). Since these excitations 
are generated from the high energy triplet components in the ground state and these 
vanish with increasing magnetic field (cf. Fig. |4.2| ), their residual spectral weight 
slowly disappears with increasing magnetization. 



5.4 Weak to strong coupling evolution 



For all the excitation spectra presented above the intrachain coupling ratio of BPCB 
7 = J||/Jx ^ 1/3.55 <C 1 was taken. For this chosen value of 7, a strong coupling 
approach gives a reasonable description of the physics. In this section we discuss 
the evolution of the spectra from weak (7 00) to strong coupling (7 0). To 



illustrate this behavior, we show in Fig. 5.8 the symmetric and antisymmetric parts of 
the correlations at = 0.25 for different coupling ratios 7 = 00, 2, 1,0. 5,0. 

At 7 ^ 00 (Fig. |5.8| l), the chains forming the ladder correspond to two decou- 
pled Heisenberg chains. In this case the symmetric and antisymmetric correlations 
are identical ^q"" = 5+" and are equivalent to the correlation 2S'+~ of the single 
chain 1 103 1 with magnetization per spin m^/2 = 0.125. A complex low energy con- 
tinuum exists with zero energy branches |[6j[2T][T03| at momenta q = ±7rm^,7r 
similar to that discussed in Sec. 



5.3.2 



In contrast, in the strong coupling limit 
(7 0) (Fig. |5.8[ 5.b) the symmetric correlations vanish and the antisymmetric 
part corresponds to the single chain correlation 2S~^~ with anisotropy A = 1/2 and 
magnetization per spin — 1/2 (see the spin chain mapping in Sec. 2.1.1| ). The 
antisymmetric part consists of a low energy continuum with branches at momenta 
q = {1 ± 2m^)7r,7r (Sec. 5.3.2). Note, that a bosonization description of the low 



energy sectors of both extreme regimes can be formulated |[3| [2H[22| (Sec. |3.2.3 ) 



^The Zeeman shift includes both the shift of the ground state (Fig. |5.7) and the shift of the excited 
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Figure 5.8: Momentum-energy dependent H correlations (S^~ (q, cj)) at = 0.25 

for different ladder couplings (7 = J\\/J±) (1) 7 ^ 00, (2) 7 = 2, (3) 7 = 1, 
(4) 7 = 0.5, (5) 7^0. The symmetric (antisymmetric) correlations with = 
(Qy — ^) ^re presented in the figures labeled by a (b). In (l,2-4.a) the vertical dashed 
lines represent the incommensurate momenta of the low energy branches of the single 
spin chain at ^ = ±7rm^ = =b7r/4 (they also correspond to the predicted momenta of 
the lowest energy excitations of the symmetric correlations | 22 |). The horizontal solid 
(dotted) horizontal lines in (2-4. a) correspond to the approximate energy J± (3J±) of the 
single triplet excitations of type (ii) (two-triplet excitations of type (iii)). The horizontal 
dashed lines in (2-4.b) correspond to the approximate energy 2J± of the excitations of 
type (ii). The dash-dotted lines in (1) correspond to the linear low energy boundaries 
of the continuum of excitations given by the LL theory applied to a single Heisenberg 
chains 7 ^ 00. The dash-dotted lines in (2-5. b) correspond to the linear low energy 
boundaries of the continuum of excitations given by the LL theory for a spin ladder with 
finite 7 (Sec.[3£3|. [Taken from Ref. j38)] 
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In the following we discuss the evolution between these two limits. In the anti- 
symmetric correlation (cf. Fig. 5.8 2-5. b) a low energy continuum exists at all cou- 
plings with a zero energy excitation branch at q = tt. These low energy excitations 
correspond mainly to the excitations with AS = AM^ = — 1. This has been pointed 
out for the weak coupling limit |[2T 1Q3| . They become the transitions \s) 
with the same quantum numbers in the decoupled bond limit. Additionally the up- 
per part of the excitation spectrum at weak coupling, which mainly corresponds to 
excitations with | |103| AS = 0, 1 and AM^ = — 1 splits from the lower part of 
the spectrum and moves to higher energy while increasing the coupling. This upper 
part evolves to a high energy excitation branch which corresponds in the decoupled 
bond limit to the \s) \t~) transition, i.e. single triplet excitations of type (ii) 
(Sec. |5. 3.3.2] approximately a|^2J^. 

The properties of the zero energy excitation branch at q = n show a smooth 
transition between the two limits |[3][2T|. For example the slope of the lower edge 
continuum which is determined by the LL velocity u decreases smoothly from its 
value for the Heisenberg chain to the lower value for the anisotropic spin chain with 
A = 1 /2 in the strong coupling limit. In contrast to this smooth change, the presence 
of a finite value of J± leads to the formation of a gap in the incommensurate low 
energy branches | 21J atq = ±7rm^ . With increasing coupling strength J± new low 
energy branches at momenta q = 7r{l± 2m^) become visible )3][22|. The weight of 
these gapless branches is very small for small coupling and increases with stronger 
coupling |3|. 

In contrast to the antisymmetric part, the symmetric part becomes gapped 
when the interladder coupling Jx is turned on. The lowest energy excitations remain 
close to the momenta q = ±7rm^ in agreement with Ref. [22) . They connect to the 
single triplet excitations of type (ii) (Sec. |5.3.3.2| ) which are approximately at an 
energy J±. While increasing 7 the higher part of the spectrum starts to separate 
from the main part and evolves to a branch of high energy two-triplet excitations 
of type (iii) (Sec. 5.3.3.3| ). These are located at approximately 3J±. Our computed 
spectra for 7 = 2, 1, 0.5 presented in Fig. |5.8[2-4.a clearly show this behavior. In 



Fig. |5.8[ 4.a the highest two-triplet excitations cannot be seen anymore since their 
spectral weight is too low. 



5.5 Influence of the weak interladder coupling on the 
excitation spectrum 

Up to now we only discussed the excitations of a single spin ladder and neglected 
the weak interladder coupling J' usually present in real compounds. 



^Note that in the strong coupling limit the energy scales set by J± and become very close, such 
that this position is in agreement with the ones previously discussed for the strong coupling limit. 
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Deep inside of the spin liquid phase, the correlations for a single ladder are dom- 
inated by high energy single or multi triplet excitations as discussed in Sec. 5.2 The 
presence of a small interladder coupling J' causes a dispersion in the interladder di- 
rection with an amplitude of order J' . This effect can be evaluated for independent 
triplet excitations using a single mode approximation However, for the com- 
pound BPCB the interladder coupling is so small that present day experiments do 
not resolve this small broadening |34, 35 1. 

In contrast in the gapless phase the effect of the interladder coupling can change 
considerably the excitations. In particular, below the transition temperature to the 
3D-ordered phase, a Bragg peak appears at = tt in the transverse dynamical func- 
tions S^^ due to the transverse antiferromagnetic long range order presented in 
Sees. |3.3| and |4.4.2] As discussed in Ref. f92l, this Bragg peak is surrounded by 
gapless Goldstone modes and it has been measured in the compound BPCB [32J 
by ND experiments described in Sec. |5.6.l] (see Fig. |5.13[ ). Additional high energy 
modes are predicted to occur in the transverse S^^ and longitudinal Sq^ \9^. It 
would be interesting to compute the excitations using random phase approximation 
analogously to Ref. | [92| in combination with the computed dynamical correlations 
for the single ladder in order to investigate the effect of a weak interladder coupling 
in more detail. However, this goes beyond the scope of the present work and will be 
left for a future study. 



5.6 Inelastic neutron scattering (INS) 

The inelastic neutron scattering (INS) technique is a direct probe for dynamical 
spin- spin correlation functions. Measurements have been performed on the com- 
pound BPCB in the spin liquid phase p4|[35| (low magnetic field) and in the gapless 
regime | 34|. Modeling the compound BPCB by two inequivalent uncoupled ladders 



oriented along the two rung vectors di^2 ( |2.14| ) (see Fig. 2.6 ) the magnetic INS cross 
section |114| is given by the formula 

^) [c(Q) • (5o+- + 5o-+) + .(Q) • + 5-+)] I (5.11) 

with c(Q) = Ei=i,2Cos2(Q • di/2) and 5(Q) = Ei=i,2sin^(Q ' di/2). Here 
Q = {Q^ Q^jQ^) = is the momentum transferred to the sample (q, are the 
ingoing, outgoing neutron momenta, respectively) and cj = E — E' is the transferred 
energy (E, E' are the ingoing, outgoing neutron energies). The correlations S^^ are 
defined at zero temperature in Eq. \5A) and evaluated at a momentum = Q • a 
along the a unit cell vector (momentum along the ladder direction) and energy uo. 
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The magnetic form factor F{Q) of the Cu^+ and the ratio q' /q are corrected in the 
experimental data. 

The INS cross section ( |5.11| ) is directly related to a combination of different cor- 
relation functions S^^ with weights depending on the transferred momentum Q and 
the magnetic field orientation. In the model definition (see Sec. |2.1| ), the magnetic 
field h is pointing along the z direction. Additionally, the two rung vectors ( |2.14| ) of 
BPCB have identical components along the unit cell vectors a and c. Hence, align- 
ing the field to the b unit cell vector and tuning Q in the a'^c'^ -plane (a"^, b"^ and c"^ 
are the reciprocal vectors of a, b and c) allows one to keep constant the prefactors in 
front of the correlations in Eq. ( |5.11| ) scanning the a"^ -momentum with the condition 
Q • di = or TT, for both i = 1, 2, to target the symmetric or antisymmetric part, 
respectively. 

We focus here on the antisymmetric part for which the low energy spectra have 
already been studied experimentally and theoretically p4| . Theoretically the focus 
so far lay on the description of the low energy excitations by the spin chain mapping. 
In contrast, we compute here the INS cross section ( |5.11| ) for the full ladder at = 
0.25, 0.5, 0.75 using the correlations presented in Sec.|5.3[ The results are shown in 



Figs. 5.9 and 5.10 and are compared to the results from the spin chain mapping. 




Figure 5.9: Theoretical momentum-energy dependent INS cross section for BPCB with 
Q • = TT (i = 1,2) and ^ = Q • a at = 0.5 in (a) a ladder system and (b) the spin 
chain mapping. The horizontal dashed lines correspond to the constant energy scans at 
uj = 0.2, 0.4 meV shown in Fig.plTl [Taken from Ref. |38 |] 



As expected from expression ( |5.11| ), the INS cross section contains the different 



excitations present in the spectra of S^^ and S^^ (cf. Figs. 5.2 b, 5.3 b and 5.4 b): 

(a) The low energy continuum originates from the transversal correlations S^^. 
It is qualitatively well described by the spin chain mapping that presents a 
symmetry with respect to half magnetization. 

(b) The continuum of excitations at energy ~ J± comes from the longitudinal 
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correlation S^^ (not present in the spin chain mapping). 

(c) The continuum of excitations at energy ~ 2 stems from the transversal 
correlation (not present in the spin chain mapping). 
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Figure 5.10: Theoretical momentum-energy dependent INS cross section for BPCB 
with Q • di = TT (i = 1, 2) and ^ = Q • a, (a) at = 0.25 (b) at = 0.25, 0.75 in 
the spin chain mapping, (c) at = 0.75. The horizontal dashed lines correspond to the 
constant energy scans at cj = 0.2, 0.4 meV plotted in Fig. 15. 121 [Taken from Ref | 38 1] 



The main features of the low energy continuum (a) are well covered by the spin 
chain mapping p4| . However, slight differences between the low energy excitations 



in the spin ladder and the spin chain are still visible in Fig. 5.10 (cf. also Sec. 5.3.2 ). 
These differences can even be distinguished in the experimental data in Figs. |5.ir 
and 5.12 where some cuts at fixed energy u = 0.2, 0.4 meV are plotted. The INS 
measured intensity is directly compared to the theoretical cross section ( |5.11| ) com- 
puted for the ladder and the spin chain mapping at = 0.24, 0.5, 0.72 convolved 
with the instrumental resolution. The amplitude is fixed by fitting one proportionality 
constant for all fields, energies, and wave vectors. 

The scans at fixed energy present peaks when the lower edge of the conti nua 
(related to the correlations S^^) is crossed (see dashed white lines in Figs. 5.9 
and |5.10[ ). As one can see, the theoretical curves for the ladder and the spin chain 
both reproduce well the main features in the experimental data and only small dif- 
ferences are present: 



- The spectral weight intensity at = 0.5 and uj = OA meV (in Fig. 5.11 b) 
is slightly overestimated by the spin chain mapping. 

- The height of the two central peaks at = 0.24 and uj = 0.2 meV (in 



Fig. 5.12 c) is underestimated by the spin chain mapping. 



Whereas the low energy excitations (a) only showed a slight asymmetry with 
respect to the magnetization, a very different behavior can be seen in the high energy 
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part (b)-(c). Indeed, the high energy part of the INS cross section (Fig. |5.10]) i s very 
asymmetric with respect to half magnetization. As we discussed in Sec. |5.3[ these 
excitations are due to the high energy triplets that can be excited in S^^ and 
(see Fig. |5.2[ b and Fig. |5.3| b) and are totally neglected in the spin chain mapping. 
Their corresponding spectral weight is of the same order than the low energy spectra, 
and thus should be accessible experimentally. It would be very interesting to have an 
experimental determination of this part of the spectrum, since as we have seen this 
part contains characteristic information on the system itself and is related to itinerant 
systems via the various mappings. 




Momentum transfer q (27ra"^) 



Figure 5.11: Inelastic neutron scattering intensity measured along a"^ of BPCB |34| with 
the momentum tt in the rung direction (Q • di = tt) at — 10.1 T (m^ ^ 0.5) and 
T — 250 mK after subtraction of the zero-field background. In each panel, fixed energy 
scans (shown by white dashed lines in Fig. |5.9| l are plotted: (a) uj = 0.2 meV, (b) uj = 
0.4 meV. The circles correspond to the experimental data. The red (black) solid lines 
are the = 0.5 theoretical data for the ladder (the spin chain mapping) convolved with 
the instrumental resolution. The shaded bands indicate the error bar in the experimental 
determination of a single proportionality constant valid for all fields, energies, and wave 
vectors. The width of these areas combines the statistics of all scans with uncertainties 
in the exact magnetization values at the chosen fields and in the convolution procedure. 
[Taken from Ref. (38J] 



5.6.1 Neutron diffraction (ND) 

The so-called neutron diffraction (ND) is an elastic neutron scattering process with- 
out energy transfer (cj = 0). Assuming a non-degenerate ground state, which is the 
case for the (weakly coupled) spin- 1/2 ladder \23) and \2A) (without frustration 
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Figure 5.12: Inelastic neutron scattering intensity measured along a"^ of BPCB f34\ with 
a TV momentum along the rung direction (Q • di = tt) at T = 250 mK after subtraction 
of the zero-field background. In each panel, cuts at fixed energy (shown by white dashed 
lines in Fig. |5.10| l are plotted: (a) u = 0.4 meV and = 0.24, (b) u = 0.4 meV 
and 0.72, (c) cj 0.2 meV and 0.24, (d) uj 0.2 meV and 0.72. 

The circles correspond to the experimental data. The solid red (black) curves are the 
theoretical data for the ladder (the spin chain mapping) convolved with the instrumental 
resolution. The shaded bands indicate the error bar in the experimental determination 
of a single proportionality constant valid for all fields, energies, and wave vectors. The 
width of these areas combines the statistics of all scans with uncertainties in the exact 
magnetization values at the chosen fields and in the convolution procedure. [Taken from 
Ref. 1^] 



(see Ref. 1 1 ])), the cross section ( [5.11| ) leads to 



dVtdE 



1 - ^ ) KQ) • + «(Q) • «(«)) 



c(Q).K^(9)r + 5(Q)-(m:^(g)r 



. (5.12) 



The magnetic orders m^^(g) = {Sl^{q)) and m^^f (q) = ^ {S^^{q))^ + {S^y{q))^ 
computed for the ground state, with momentum q (qy) along the leg (rung), are par- 
allel and perpendicular to the magnetic field orientation, respectively. Hence, fol- 
lowing Eq. ( |5.12| >, at low temperature, ND provides a quantitative measurement of 
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the order parameters. It has been used to measure on BPCB the longitudinal mag- 



netization (Fig. 4.9 ) and the transverse long range order in the 3D phase. For the 
latter, Fig. |5.13 a shows the intensity of the ND experiments in the gapless regime 
versus the momentum Q along a"^. The appearance of a Bragg peak at the momen- 
tum g = Q • a = TT for temperatures below Tc{h^) is the signature of the existence 
of a staggered long range order. Following the ND cross section ( |5.12| ), the inten- 
sity of the peak is proportional to (mj^(7r))^ and is used to extract the 3D staggered 
transverse order m^{h^) shown in Fig. |4.12[ In Fig. 5.13 b, the intensity of the peak 
plotted versus the temperature at fixed magnetic field clearly shows the onset of 
the 3D order at Tc{h^). These ND measurements from Ref. |32| are used to extract 
the magnetic field dependence of the critical temperature Tc{h^) shown in Figs. 2.4 
andl4l2l 
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Figure 5.13: Neutron diffraction measurements on BPCB at several temperatures T and 
applied magnetic fields B (h^ in our notation) in the gapless regime, (a) Momentum 
scan Q along a"^ of BPCB across an antiferromagnetic Bragg peak after subtraction of 
a flat back-ground measured in the spin liquid phase aih^ = 6 T and T = 63 mK. (b) 
Temperature dependence of the Bragg peak intensity, demonstrating the onset of 3D long 
range order at Tc{h^). SoHd lines are fits using the 3D-XY exponent 2/3 ^ 0.7 |[4l |115| . 
[Taken from Ref. f32l] 
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CHAPTER 6 



Conclusions and perspectives 



6.1 Conclusions 

In this work we have investigated the thermodynamic and dynamic properties of 
weakly coupled spin ladders in a magnetic field. Combining a LL analytical theory 
and DMRG numerical techniques, both generalized using a mean field approxima- 
tion while considering the interladder coupling, we were able to explore characteris- 
tic physical properties of the system in all the regimes of the phase diagram shown 
in Fig. |2.2| In addition to the theoretical analysis we compared our findings to exper- 
imental results on the compound BPCB ((C5Hi2N)2CuBr4) and obtained excellent 
agreement. This compound appears to be a very good realization of such ladder 
systems in which other possible effects such as frustration, anisotropic interaction 
or Dzyloshinskii-Moriya term are very small. The excellent low-dimensionality of 
this compound and the ideal range of its coupling parameters lead to a clear sep- 
aration of the energy scales between the different phases which were all accessed 
experimentally. 

As thermodynamic quantities we computed the magnetization and specific heat 
of the system as a function of temperature and magnetic field. The extension of 
the DMRG technique to finite temperature allowed us to compute these quantities 
with an excellent accuracy. In the gapless phase the low energy part of the specific 
heat agreed well with the prediction of the LL theory. At higher temperatures, the 
numerical solution was needed to capture the precise structure of the peaks in the 
specific heat, that reflect the presence of the excited states in the ladder. The com- 
parison of the theoretical calculations with the measured magnetization and specific 
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heat proved to be remarkable. This good agreement confirms that the ladder model 
is indeed a faithful description of the compound BPCB. It also gives direct access, 
via the extrema in magnetization and peaks in the specific heat to the approximate 
region of applicability of the Luttinger liquid description. For BPCB we found an 
extended region lying approximatively one order of magnitude above the 3D order 
transition temperature. The large domain of validity of the LL description for BPCB 
opened the way for a more detailed experimental investigation of this regime. 

For the low energy dynamics we used a combination of the numerical techniques 
to determine the Luttinger liquid parameters and then the analytical description based 
on Luttinger liquids to compute the dynamical spin-spin correlation functions. This 
allowed us to extract the NMR relaxation rate and the one dimensional antiferro- 
magnetic transverse susceptibility. If the ladders are very weakly coupled, which is 
the case in the considered material, the divergence of the susceptibility leads to a 
three dimensional antiferromagnetic order at low temperatures in the direction trans- 
verse to the applied magnetic field. We computed this transition temperature and the 
order parameter at zero temperature. Comparison with the measured experimental 
quantities both by NMR and neutron diffraction proved again to be remarkable. This 
excellent agreement between theory and experiment for these quantities as a function 
of the magnetic field allows to quantitatively test the Luttinger liquid theory. It shows 
that several different correlations are indeed fully described by the knowledge of the 
two Luttinger liquid parameters (and the amplitudes relating the microscopic opera- 
tors to the field theoretical ones). This is something that had not really been tested 
previously since either the microscopic interactions were not known in detail leaving 
the Luttinger parameters as adjustable parameters, or only one correlation function 
could be measured in a given experiment, not allowing to test for universality of the 
description. 

We also gave a detailed analysis of the dynamical spin-spin correlations, for T = 
using the time dependent DMRG method, for a wide range of energies and all 
momenta. The excitations reveal a lot of important information on the system and 
are well suited to characterize it. In particular we showed the interesting evolution 
of the excitations in the system with the magnetic field and the coupling strengths. 
Quite interestingly the intermediate energy part can be related to the excitations of a 
t-J model and shows thereby features of itinerant systems. We also showed that the 
dynamical correlations of the ladder posses characteristic high energy features that 
are clearly distinct from the corresponding spectrum for spin chains. 

The numerical calculation is efficient for the high and intermediate energy part 
of the spectrum for which the Luttinger liquid description cannot be applied. We 
showed that the two methods, numerics and LL have enough overlap, given the ac- 
curacy of our calculation so that we can have a full description of the dynamical 
properties at all energies. This allowed us to use each of the method in the regime 
where it is efficient. In particular, in this thesis we did not push the numerical cal- 
culations to try to obtain the exact behavior at low energies or the very fine structure 
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effects, but focused on the high and intermediate energy regime suitable for the ex- 
isting and future INS experiments. We used the analytical description coupled to 
the numerical determination of the Luttinger parameters to obtain an accurate low 
energy description. We made the connection between our results and several analyt- 
ical predictions. In particular at intermediate-low energy our calculation agreed with 
the Luttinger liquid prediction of incommensurate points and behavior (divergence, 
convergence) of the correlations. 

We compared our numerical results with existing INS data on the compound 
BPCB and found excellent agreement. It is rewarding to note that the resolution of 
our dynamical calculation is, in energy and momentum, at the moment better than 
the one of the experiment. The comparison between theory and experiment is thus 
essentially free of numerical errors and thus allows for precise tests of the models. 
Given the current resolution of the INS experiment it is difficult to distinguish in the 
low energy part of the spectrum the difference between the dynamical correlation of 
the true ladder and the one of an anisotropic spin 1/2 system, which corresponds 
to the strong rung exchange limit. Additionnal experiments would be desirable in 
this respect. An alternative route is to probe experimentally the high energy part of 
the spectrum, since the predicted high energy excitations contain many characteristic 
features of the underlying model. 

In connection with the compound BPCB and also on the conceptual side, there 
are several points which remain to be investigated: 

• An improvement of the description of the quasi one dimensional systems, by 
including in a mean field way the effect of the other ladders in the numerical 
study, has been partly performed in this thesis. Nevertheless the extension to 
the dynamical quantities remains to be done. This is specially important close 
to the quantum critical points hd and hc2 where the interladder coupling be- 
comes crucial and the system undergoes a dimensional crossover between a 
one dimensional and a higher dimensional (three dimensional typically) be- 
havior. Understanding such a crossover is a particularly challenging question 
since the system goes from a description for which a picture of essentially free 
fermions applies (in the one dimensional regime) to one for which a descrip- 
tion in terms of essentially free bosons (the three dimensional regime) applies. 



As discussed in Sec. |4.5| small discrepancies between the 3D order parameter 
measured on BPCB at finite temperature and the computed mean field value 
at zero temperature remain unexplained. In order to clarify these deviations 
an investigation of the temperature effects on the order parameter would be 
necessary. In addition a more precise description of the 3D structure going 
beyond the mean field approximation, already partially done in Ref. 1 32 1 using 
a quantum Monte-Carlo technique, including eventually frustrating interladder 
couplings would probably help to understand the discrepancies. 



As briefly mentioned in Secs.|2.2|and 4.5 recent ESR measurements on BPCB 
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pQ| confirmed by a refined theoretical analysis |116| have pointed out small 
anisotropics in the exchange couplings of this compound. These deviations 
are very small and thus do not change in an essential way the results for 
the thermodynamic quantities and the dynamical ones presented in this thesis. 
However at lower energy scales and in some other observable they will lead 
to deviations compared to the ideal Heisenberg model which deserve further 
investigations. These investigations could also help to understand the small 
deviations of the experimental measurements on BPCB from the theoretical 
predictions discussed in Sec. |4.5| 

An extension of the dynamical results to finite temperature would be desirable. 
This could be used to study different effects such as the interesting shifts and 
damping of the triplet modes with temperature that have been observed in three 
dimensional gapped system [ 1T7| . Similar effects have been also predicted in 
the excitation continuum of spin- 1/2 chains | |78J . 

Recently, the presence of non LL edge singularities due to the band curvature 
of the dispersion relation have been pointed out in spin- 1/2 chains |106 - 108| . 



A very precise numerical analysis of the dynamical correlations close to the 
continuum threshold of the spin- 1/2 ladder excitations would be required in 
order to investigate such fine deviations from the LL description. 



6.2 Outlook 

Our combined numerical-analytical methods approach is quite general and could 
be extended to explore many other systems. In effect, various ladder structures in- 
cluding frustration, dimerization or long range interaction remain to be explored. 
Focussing on the dynamical properties, which are difficult to access, our theoretical 
approach could be extended to study the excitations of these more complex systems. 

Motivated by the impressive description of BPCB with a simple spin- 1/2 ladder 
model several other ladder compounds with different coupling ratios 7 have been 
synthesized. One of the most promising compound is (C7HioN)2CuBr4 (DIMPY) 
|118| . This material has been investigated experimentally using INS and specific 
heat measurements 1 119[|12Q| . These first experiments agree on the ladder structure. 
Interestingly, DIMPY has a coupling ratio 7 ?^ 2 and thus belongs to the opposite 
weak coupling limit compared to BPCB. A more detailed theoretical analysis of this 
weak coupling regime would thus help to characterize this compound. 

An other possible direction would consist in looking at the effects of disorder or 
doping on the ladder compounds. The former could lead to interesting physics such 
as a Bose glass phase. Indication of such a phase has been already observed exper- 
imentally in other magnetic compounds 1 121 1 and investigated theoretically |^122j- 
|125| . In contrary, the latter would in principle lead to superconductivity |126| . Fur- 
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ther theoretical investigations in connection with these two experimental realizations 
are of course strongly required for a correct interpretation of the experiments. 
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APPENDIX A 



Strong coupling expansion of a single spin ladder 



In this appendix we show how the spin ladder Hamiltonian ( |2.3| ) at strong coupling 
(7 <C 1) can be expressed in bosonic operators acting on single bonds introduced in 
Ref. 1 127 1 . This representation classifies the excitations with respect to their energy. 
We first derive perturbatively an effective system based on this Hilbert space orga- 
nization by energy sectors. We introduce the Schrieffer- Wolff transformation that 
maps the physical system to the effective one, and approximate the effective system 
using a strong coupling expansion. We evaluate the rung densities of the ground state 
in the spin liquid, and derive an effective theory for the gapless regime. Furthermore 
we evaluate the corrections of the LL parameters from the spin chain mapping. 



A.l Strong coupling expansion 

The four-dimensional Hilbert space on each rung / is spanned by the states \s), 

and (cf. Eqs. ( [2.7| ) and ( [2.8| )), obtained by applying the boson creation 
operators sj, tj tj and tj_ to a vacuum state. A hardcore boson constraint 
applies on each rung i.e. 

Qi,s + Ql,+ + Qi,o + Qi- = 1 (A.l) 

where gi^s = sjsi and gi^k = tj ^ti^k are the density operators and A: = ±, 0. 

While the Hamiltonian on the rung H±^ \2A\ is quadratic in the boson operators 

L 

H± = Y1 - hyj±)Qi,+ + Qi,o + (1 + hyj^)gi^_] - -L, (A.2) 
1=1 
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the chain Hamiltonian H\\ \2A\ is quartic, and its structure is quite complex. The 
advantage of the boson representation reveals itself when considering the case of 
small 7. In that case we perform a Schrieffer- Wolff transformation of the spin ladder 
Hamiltonian \23) 

i^eff = e'^^i^e"'^^. (A.3) 
The Hermitian operator A can be expanded in powers of 7 

A = Ai + 7A2 + • • • . (A.4) 

Thus i^eff can be written in orders of 7 as 

Jj^i/eff = H^^ jH^^^ + 72/7(2) ^ . . . ^ 

where 

H^'^ =H\\^i[Ai,H^], (A.6) 
H^^^ =i[A2,H^] - (A.7) 

etc. Through this expansion, the unitary transformation e*^^ can be perturbatively 
determined computing the Af^ recursively in order to eliminate the transitions be- 
tween the energy sectors of excitations in i^eff. Since the first term J±Hj_ in Eq. (|A.5|) 



leads to a separation of excitations on the order of the energy scale J± (cf. Fig. 2.2 a) 
the decoupled bond limit provides an effective Hilbert space that contributes to each 
energy sector. The second term Jyi^^^^ causes broadening of these bands on the 
order of Jy and can induce a complex structure within the energy bands. To obtain 
the desired expansion up to the first order in 7 we choose 

^1 = ^ ^ 44+1 {ti,oti+i,o - ti,+ti^i- - ti-ti^i^^) + h.c, (A.8) 
I 

where h.c. stands for the Hermitian conjugate. 

A.2 Rung state density in the spin liquid 

In the spin liquid phase, the decoupled bond limit provides the effective ground state 
|Oeff) = 1^ • • • 5) which is related to the physical ground state by 

|0)=e-^^^|Oeff). (A.9) 

So the triplet density of the ground state (pk) (with A: = ±, 0) is given by 

(Pk) = {Qi,k) = {s--- s\e'^^gi^ke-'^^\s • • • 5). (A.IO) 
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Using Eq. ( |A.8| ), and keeping only the non- vanishing corresponding terms in ( |A.10| ) 
up to second order we get 



(A.ll) 



In the case of the compound BPCB (see Sec. 2.2) this expansion gives {pk) = 0.01, 
and due to the hardcore boson constraint (Eq. ( |A.l| )) (ps) = {qi,s) — 0.97. Even 
though we took into account only the first order term for A in Eq. ( |A.4| ), this ap- 
proximation of the triplet density differs from the direct numerical computations (in 
Fig.|42])by only - 20%. 



A.3 Effective Hamiltonian in the gapless regime 

The first order term H^^"^ of the effective Hamiltonian ( |A.5| ) is computed substitut- 
ing ( [A.8| ) into ( |A.6| ). This leads to H'^^^ in the form 

k=0 I 

^ V ' 

where 

^'oj = + ^Qw,+Qi,+ + h.c, (A.13) 

^ij = + 4+i,+4,o^^+i,o^^+ + h.c, (A.14) 



+ 4+i,o4,o^^+i -^^+ + 4+i,o4,o^^+i,+^^- + h.c, (A.15) 

^3} = 4+i,o4,-^^+i,-^^o + h.c, (A.16) 

and 

n^li = Qi+i,-Qi,- (A.17) 

Here we regrouped the terms such that each J^H^^ acts on the corresponding en- 
ergy sector kJ^, k = 0, 1, . . . , 4 in the gapless regime. Note that in each sector 
Ai = such that to the given order in 7 the Hamiltonian ( |2.3| ) corresponds to the 
effective Hamiltonian. 



85 



A. STRONG COUPLING EXPANSION OF A SINGLE SPIN LADDER 



A.3.1 Low energy sector 

When focusing on the low energy sector, the \s) and the modes dominate the 
behavior and we can assume a vanishing density of |to) and triplets. Thus the 
hardcore boson constraint ( |A.1| ) simplifies to 

Qi,s^Qi,+ = l (A.18) 

and the rung Hamiltonian ( |A.2| ) to 

H^ = {1- h'/J^) Qi^^ - -^L. (A.19) 

I 

Furthermore the only contribution to the first order term in 7 comes from 
Taking this into account we obtain from Eq. \A.5) for the Hamiltonian \23\ in the 
lowest energy sector 

H = J^H^^J\\H^^\ (A.20) 

where is given by Eq. ( |A.19| ) and H^^ by Eq. ( |A.13| ). Following Ref. |^|, we 
map the two low energy modes onto the two states of a pseudo spin (Eq. \2.9\ ) and 
replace the boson operators 5^ and t^_^ with the spin- 1/2 operators 

St=ti+su Si=s\ti,+, St = Qi,+ -\. (A.21) 
The effective Hamiltonian is the XXZ spin-1 /2 chain Hamiltonian, Eq. ( |2.11| ). 

A.3.2 Sector of energy 

The effective Hilbert space of the J±_ energy sector corresponds to a single \ t^) triplet 
excitation lying in a sea of singlets \s) and triplets The effective Hamiltonian 
up to first order in 7 is given by 

H = J^H^ + J|| (i/^'^ + h[^^^ . (A.22) 

The excitation \t^) can be interpreted as a single hole excitation in a spin chain 
formed by \s) and Each rung state of this sector is identified with 

\l) = \s), n) = \t+), \0) = \t°). (A.23) 

In this picture the Hamiltonian ( |A.22| ) can be mapped onto the anisotropic t-J 
model 

H,.j = Hxxz + i^t + i^s-h + e. (A.24) 
where e = ( + /i^)/2 is an energy shift. 
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The hopping term 

l,a 

Stems from the term J\\H[^-' in Eq. ( |A.22| ). Here cj ^ (ci^a) is the creation (annihila- 
tion) operator of a fermion with pseudo spin a = f, | at the site I. Note that although 
we are dealing here with spin states, it is possible to faithfully represent the three 
states of each site's Hilbert space (\s), using a fermionic representation. 

Additionally, a nearest neighbor spin dependent density-density term arises 



Jii V r 

I 



(A.26) 



Here ni^h is the density operator of the hole at site I. This term stems from the 
nearest-neighbour interaction between the triplets, i.e. the second term in Eq. 
( |A.13| ). Mapping this term onto a spin chain in the presence of a hole leads to an 
interaction term between the hole and the spin up state. Note that this is in contrast 
to the usual mapping onto a spin chain without a hole in which case the term would 
only cause a shift in energy. 



A.4 Second order perturbation and Luttinger liquid 
parameters 



The second order term H'^'^^ ( |A.7| ) of the expansion ( [A.5| ), contains a huge amount 
of terms. Nevertheless considering the low energy sector, only the following terms 

^0^^ = - ^ X] Qi,sQi+i,s + h.c.^ (A.27) 



are important. The first term in Eq. ( |A.27| ) is a singlet density-density interaction 
which can be absorbed in to th e coupling of the XXZ chain, and the second term 
is a conditional hopping In order to study the effects of i^Q^^ on the LL 



parameters u and K (see Fig. 3.4), we first replace the boson operators with the 
spin- 1/2 operators (Eq. ( |A.21| )). So the Hamiltonian ( |2.3| ) in the low energy sector 
becomes 



H = H^ 



XXZ 



-E 



const, (A.28) 



where i^xxz is the XXZ spin- 1/2 chain Hamiltonian (Eq. ( |2.11| )) with the corrected 
parameters 

A(2) = ^-^7 , h^''^^ =h^-\j..l^ (A.29) 
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up to the second order in 7. For the BPCB parameters (Eq. 2.13 ) A*^^^ = 0.4 instead 
of A = 0.5 for the spin chain mapping (first order approximation). The LL parame- 
ters K and of i^xxz with the anisotropy A*^^^ are computed, and we treat the 
conditional hopping term by approximating 1/2 — Sf = 1 /2 — (mean field ap- 
proximation). The remaining term is then bosonized using the expression ( |3.26| ) for 
S^{r = I). It leads to the corrected LL parameter u and K of the Hamiltonian ( |A.28| ) 
through the relations 



uK = uK ^27r j^J^Aa: {1/2- 
u/k = u/K 



m 



(A.30) 



The corrected u and K are plotted in Fig. 3.4 and clearly show the asymmetric signa- 



ture of the full ladder parameters induced by the conditional hopping term in ( |A.28| ). 
Note, that the lack of convergence K ^ 1 when = rfi^ + 1 /2 ^ is obviously 
an artifact of the mean field approximation 1/2 — S[ = 1 /2 — . 
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APPENDIX B 



Static correlations in a finite size Luttinger liquid 



In this appendix we recall the analytical expressions of the static correlation func- 
tions of the spin- 1/2 chain and the spin- 1/2 ladder p.3| ) models computed for 
finite size systems using the LL theory (Sec. \3.2\ in Ref. |28|. These correlations 
allow for a more precise determination of the LL parameters by comparison with the 
numerical DMRG computations (see Sec. |3.2.2| ) taking into account the boundary 
effects neglected in the infinite size expressions (Eqs. |3.3l] and |3.32"| ) but which are 
obviously present in the numerical computations. 



B.l Static correlation functions 

In the following, we give the relations used for the spin chains only, since from these 
the relations for the spin ladders can be easily inferred using the spin chain mapping 
(Sec. [2TT] ) 

{S^^^Sl^)^2{SfSf.), (^f,o^f^o)^(^f^^) + ^((^^^ (B.l) 

The other correlations {Sf^Sf, q) and {S[^S[, ^) are not present in the spin chain 
approximation. Using a weak coupling approach | [22| or a direct DMRG computa- 
tion j28| these correlations are shown to vanish exponentially when |/ — oo 
leading to the gapped spectra of the corresponding energy-momentum correlation 



functions shown in Figs. 5.3 a, 5.4 a and 5.2 b. 
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B. STATIC CORRELATIONS IN A FINITE SIZE LUTTINGER LIQUID 



The LL analytical expressions of the correlation functions (SfSf,), (SfSf,) and 
the local magnetization (Sf) for a system of length L are 



/i/4if(20/i/4if(2r) 



^x(-l) 



l-l' 



h/2K{l-l')h/2K{l + l') 

(-l)'cos(qr) {-if cos{ql) 



^2A,B, sgn{l - I') 

' fK{2i)fKm 



fKm 



r f, , f2K{l - I') 

cos[q{l + I )\ 



(B.2) 



27r 



(-l)'sin(gO {-if sm{ql') 



27r [27r ' ^ \^ /k(20 
K f I 1 

2^ KJaT^^WTV) 



fK{2l)fK{2l') 

KV2A~^ f{-iycos{ql) 
/if (20 



cos[q{l - I )] - cos[q{l + I )] 1 



'f2K{i-n 

[g{l + l')+g{l-l')] 



i-iy' cos(gl') 



/k(2?') 



[ff(l + n-<7a-0] 



(5f ) = 



27r 



■(-l)'sin(90 
/if (20 



(B.3) 
(B.4) 



with 

/a (a;) 



2(i + l) . 

^ ^ sm 



2{L+1)J_ 



2{L + 1) 



cot 



2(L+1)^ 

(B.5) 

and q = 2-KLrh^/{L + 1). We can verify that in the limit of infinite system size 
(L —■ oo) and far from the boundaries (\l - L/2\ < L and \l' - L/2\ <C L) the 
LL correlations Eqs. ( |B.2| i and ( |B.3| ) simpUfy to the well known power law decay 
Eqs. ( |3.31| l and ( |3.32[ >, and the local magnetization ( |B.4| > becomes constant, (S'f ) = 
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